Rac'ti PCT/PTO ! 7 DEC 2004 

WO 03/107137 PCT/US03/19179 



TITI/E OP INVENTION: METHODS, SYSTEMS AND COMPUTER PROGRAM PRODUCTS TO 
FACILITATE THE FORMATION AND TRADING OF DERIVATIVES CONTRACTS 



BACKGROUND OF THE INVENTION 

Financial derivatives, also known as contingent claims, arc special types of contracts used by financial risk managers 
to hedge against the fluctuations of more basic undcrlyings throughout time. Throughout this document, the term 
underlying is defined as in FAS Statement 133, appendix A, paragraph 57, a. In the recent years, derivatives have 
become increasingly important in the field of finance. According to the latest report from the Bank for Internationals 
Settlements (BIS) on Derivatives activity, the notional amount on all derivatives positions held by trading institutions 
stood at over 187 trillion dollars at the end of year 2000 (this value represents the sum of contracts traded through 
exchanges and over the counter). The financial industry has responded to this growth in the derivatives market by devel- 
oping new methods and systems to price and hedge derivatives contracts and to facilitate trading in derivatives contracts. 

The following problems arc known to limit the accuracy in present approaches to price and hedge derivatives contracts, 
hence the ability to most efficiently and effectively trade derivatives contracts: 



a) a) The current approaches to valuing derivatives make the assumption that hedging is done in continuous time, 
while in practice, participants hedge their positions in discrete time increments. Therefore, the existing approaches 
used in the financial industry or suggested by the academic community do not accurately account for this inherent 
discrepancy and arc limited in accuracy when applied to real markets and the pricing of derivatives. 

b) Current approaches to valuing derivatives assume that markets arc frictionlcss — bid/offer spreads arc reduced to 
zero, order size docs not affect price or inventory, slippage effects and credit risk arc non-existent. In reality, however, 
markets arc not frictionlcss. While theoretical methods have been suggested to address this problem, they still 
inadequately reflect real market situations. This forces practitioners to use sub-optimally efficient methods. 

c) The current methods for derivatives pricing and hedging depend on the dynamic of undcrlyings model, thus creating 
a model risk that can be very costly if not accounted for in practice. 

d) The current methods for derivatives hedging use parameters known as Greeks. Greeks arc obtained by differentiating 
price with respect to various other model parameters. The implicit assumption behind using Greeks is that derivatives 
prices arc polynomial functions of these other model parameters. This is not true and therefore the use of Greeks in 
hedging creates an additional source of approximation error that can be costly in practice. 

c) The current methods for derivatives pricing and hedging assume that the undcrlyings move in infinitely small incre- 
ments while in practice for all markets there is a minimum increment size often referred to as a tick or a pip. 

f ) When not assuming continuous time, the current methods for derivatives pricing deal only with a single future period 
or single underlying scenario. In reality, there arc multiple periods and undcrlyings to consider. 

g) The existing methods that do attempt to address all these shortcomings arc intractable in practice and thus fail to 
provide the benefits they were set up to yield. 

A vivid example of the potential cost and risk to the financial system of the inaccuracies discussed above is the collapse 
of the hedge fund Long Term Capital Management (LTCM), where the Nobel Prize winning fund managers relied on 
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a model from which the market deviated. The real market deviation from the model led to an exposure of about one 
trillion dollars, prompting the federal reserve to engineer a bailout to avert a failure that would have otherwise disrupted 
the whole U.S. financial system and could have easily have extended to all international markets (Sec for example "When 
Genius Failed: The Rise and Fall of Long-Term Capital Management" by Roger Lowcnstcin ISBN 0-375-75825-9 [46]) 

There is a substantial body of academic and published patents addressing derivatives trading issues. The relevant 
patents and patent applications found can be classified into five categories, as described below. 

1. Patents that provide methods and systems to price specific types of derivatives* In this category, we can 
note U.S. 4,642,768, U.S. 5,799,287, JP 2001067409 

2. Patents that provide methods and systems to automate the derivatives pricing process. In this category, 
we can note U.S. 6,173,276, U.S. 5,692,233, and patent applications US20020010667 and US20020103738 

3. Patents that provide methods and systems to speed up the derivatives pricing process. In this category, 
we can note U.S. 5,940,810 and U.S. 6,058,377. 

4. Patents that provide methods and systems to better hedge derivatives or manage the risk of derivatives 
books. In this category, we can note, U.S. 5,819,237, US6,122,623 and patent applications US20020065755 and 
WO0133486 

5. Patents that provide methods and systems to more efficiently trade specific types of derivative con- 
tracts. In this category, we can note U.S. 4,903,201, U.S. 5,970,479 , U.S. 6,421,653, U.S. 6,317,727 and U.S. 
6,347,307 

All those patents while presenting benefits over their prior art face limitations addressed in this invention. 

In a recent article in the Journal of Finance on patenting in the field of finance methods and formulas, the author 1 
found problematic the failure to cite academic research in prior art reviews of patent applications or granted patent. 
Our detailed description thus start by reviewing reviewing the academic literature. 



1 Josh Lerner n Where does State Street Lead? A First Look at Finance Patents, 1971-2000" THE JOURNAL OF FINANCE VOL LVII, 
NO 2 APRIL 2002 [45] 
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BRIEF SUMMARY OF THE INVENTION 

This invention relates to methods, systems and computer programs to facilitate the formation trading and risk manage- 
ment of derivatives contracts on one or more underlying. 

It provides a new framework and significant improvements in the following areas: 

1. The decomposition of any derivatives contract into fundamental building block structures called basis instruments 
that arc financial contracts in a multi- period, multi-securities market. 

2. The incorporation of supply and demand price sensitivities in the pricing of derivatives contracts. 

3. The incorporation of credit risk in the pricing of derivatives contracts. 

4. The development of methods, systems or computer program products for the accounting of derivatives contracts 
in compliance with FAS 133 or IAS 39. 

5. The development of methods, systems or computer program products for the pricing of derivatives contracts. 

6. The development of methods, systems or computer program products for the risk management of derivatives 
contracts. 

7. The development of methods, systems or computer program products for the trading of derivatives contracts 
whether in organized exchanges or in over-the-counter (OTC) markets. 
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B RIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS 

Fig 1 gives an overview of the preferred embodiments of the invention. 

Fig 2 describes the sequence of contractual agreements and cash flows in a BIC. 

Fig 3 describes the format description of derivatives contract in the DCWBSOF format and how it is compressed in the 
DCWOF format. 

Fig 4 describes the steps of the iterative decomposition process, how derivatives contracts and BICs arc itcrativcly 
combined to yield the price of the derivative contract. 

Fig 5 describes an embodiment of the present invention in a pricing system. 

Fig 6 shows a scrcenshot of an embodiment of the present invention implemented in an online pricing system wherein 
any user specify and name their derivatives contract payout payment function (paoff function) in a functional format as 
well as their BICs in a functional format and system parameters such as time and space increments. 

Fig 7 shows how BICs functional format is input in an HTML page. 

Fig 8 shows how the derivatives contract in DCWBSOF or DCWOF is input in an HTML page. 
Fig 9, 10,11,12 & 13 show an embodiment of the present invention in an exchange system. 

The exchange trading system of Fig 9. comprises three types of players: 

• The price takers operating 1 

• The market makers operating 3 and 4 

• The exchange system management operating 2 and responsible for control and compliance operations, clearing 
and settlement verifications and more globally exchange system risk management. 

1 contains an input interface 11 through which the price taker may describe a product of interest for pricing informa- 
tion 15 or put an order 16 on an identified product. That input interface 11 may be a keyboard, a mouse, a pad, a 
microphone or any sensor capable of understanding messages transmitted by human senses; 1 also contains an output 
interface 12 that may be a screen, a speaker, a printer, or any device transmitting signals decipherable by human senses 
and understandable by the human brain. 1 also contains an authentication process 13 that uses various authentication 
algorithms to coiifirm the identity of the price taker prior to authorizing any operation with the exchange system. The 
pricing information 15 as well as the order confirmation 17 or the trade confirmation 18 arc channeled through that 
output interface 12 to the price taker senses. Through the input interface 11, the price taker sends information relating 
to the product of interest as in 19 of Fig. 1. As further specified in Fig. 2, the product of interest is specified in the 
form of a function f as / (5j , • - • , ST\ • • • , S* • • • , S£) . The function / may be described in a computing language that 
could be easily understood by price takers and translated by the server of the exchange system in 28 to be ready for 
decomposition in basis instruments. 



WO 03/107137 



5 



PCT/US03/19179 



2 contains compliance and control systems 27 that check and authorize operations performed by the price takers or the 
market makers. As further described in Fig 3., 28 performs, if necessary, the transformation of relation (26) described in 
our detailed description that takes as input a function of unspecified parameters /3» and H» as well as the S{ *s and returns 
as output to 21, a function of solely thoS^'s. 21 further taking as input the most competitive basis instruments prices 
with the references of the market makers quoting them as provided by 22, uses the iterative decomposition algorithm of 
relations (20) to compute the price of the desired product which is then returned to 15. The basis instruments prices in 
22 arc taken from individual market makers quoted basis instruments prices in 25 and 26 for example and by selecting 
the most competitive of each of the basis instruments with the references of the quoting market maker. If the price 
taker now decides to put an order on the product, this indication is channeled through 16 to 21 which then sends the 
order confirmation through 17 to the price taker, and the references of the price taker as well as the market makers arc 
sent to 24 for clearing and settlement. Trade Confirmations arc then sent through 18 to the price taker and through 32 
to each relevant market maker with information on the nature of the basis instrument contracted. 

The exchange system may operate as the counterparty in each transaction cither with the market makers or the price 
takers. In this instance, the exchange system would act as the guarantor of each obligation cither to the market makers 
or the price takers. It would then put in place a credit management function to hedge against the risk of default of any 
counterparty. This could be done through traditional means such as margin requirements or through credit derivatives 
by allowing basis instruments based on undcrlyings described in our detailed description to trade and by entering into 
contracts in which the credit risk is shifted to the most credit worthy players. The exchange system may also decide 
to let each party bear the credit risk of his or her counterparty and act solely as means of exchange, but then credit 
exposure management would be a very complex task for each and every counterparty to handle. If an exchange system 
is setup to be the counterparty in each transaction and if the credit risk function is managed properly, with the exchange 
set to have the highest credit rating possible so as to be considered the risk-free counterparty of reference, this would 
give such an organized exchange a competitive advantage over OTC derivatives markets. If furthermore it decides to 
hedge its exposure by entering into credit derivatives contracts, this would create new markets for credit derivatives. 

Another critical feature of the exchange system for computational purposes is that when the number of possible states 
of the selected random variables grows large and the number future periods grows as well, tho computational cost to 
price a derivatives contract becomes at least of the order of 0(s nm ). So an actual system would include some dimension 
reduction variables transformations to reduce the number of variables on which the various derivatives contracts prices as 
well as the various basis instruments price depend. Also, in order to reduce the number of states, some actual algorithms 
would sample only some limited number of states to infer the price of larger combinations of basis instruments prices. 

3 and 4, which arc just 2 examples of systems operated by a potentially larger number of market makers, contain an 
input interface 35 as in 11 and an output interface 36 as in 12. It also contains an authentication process 37 for market 
makers as in 13 for price takers. Bach market maker selects a model and model parameters 34 that can be processed in 
31 to yield the market maker's quoted basis instrument prices that arc subsequently sent to 25. 34 is further detailed in 
Fig. 3 and shows that the model and model parameters description function is divided into a counterparty credit risk 
model 341, models of the stochastic process followed by the undcrlyings 342, models of the scaling density functions 
343, and models for delimiting the range of basis instruments quoted 344 , all this combined as further defined in our 
detailed description to yield in 31 the quoted basis instruments prices. 3 and 4 also contain in their Assets and Trading 
Information Database 33, an inventory of the market maker's position in all basis instruments held. Trades Confirmation 
and settlement information from 24 arc sent to the market makers through 32 where the information formatted for the 
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output interface. All relevant information storing is done in 33. 

While the drawings of figure 9 only involve a single price taker, only two market makers and only a single Exchange 
System Management administrator, it is fairly obvious this is to case presentation and that an actual such system may 
involve a larger number of price takers, market makers and exchange system management personnel. 

Fig 14. describes an embodiment of the present invention in a risk management system. 

The description of the risk management system is based on our decomposition of any derivative in basis instruments and 
our method for obtaining the implied conditional probability density functions from market basis instruments prices. 
The Risk Manager communicates with the Risk Management or Decision System 6 via an input interface 68 as 51in 
5, 11 in 1 or 35 in 3 or 4, and an output system 67 as 52 in 5, 12 in 1 or 36 in 3 or 4. The Authentication Process 
in 65 is as 53 in 5, 13 in 1 or 37 in 3 or 4. The Risk Manager sends portfolio position requests to 66 which process 
them in 64, which then send a request to 63 to check the inventory of basis instruments positions held, another request 
may also be sent to 61 to obtain the market basis instruments prices. Using our formula [16] obtained in our detailed 
description to extend the Brccdcn Litzcnbergcr formula, wc obtain the conditional probability density function in 62, 
which together with the basis instruments help obtain the value of the portfolio held, compute various pro-forma or 
non pro-forma reports which help better apprehend the risk and opportunities of the position held. For example, value 
at risk scenarios can be generated with exactly market-accurate probability of occurrences directly linked to the actual 
costs of available hedges. 

Fig 15 & Fig 16 show how an embodiment of the present invention enables Value At Risk (VAR) risk management. 

Fig 17 shows how a portfolio managed in a VAR framework potentially exposes to potentially unbounded losscs(yl) in 
likely scenarios ( state between xll and x!2) that may fall out of the confidence interval. On the contrary, an embodiment 
of our invention enables risk management in which a maximum worst case scenario loss can be ascertained with certainty 
by enabling precisely designed contracts to achieve that goal. 

Fig 18 shows how computer implemented methods may be used to reduce computational time in the specification of 
BICs and the resulting derivatives price computation. 

Fig 19 shows how computer implemented methods may be used to reduce computational time in the specification of 
derivative contracts and the resulting derivatives price computation. 
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DETAILED DESCRIPTION OF THE INVENTION 

1 Introduction 

The present invention introduces the notion of basis instruments as the building blocs for pricing or hedging any deriva- 
tive security d efin ed in the most general sense. It describes a decomposition formula that precisely shows how any 
derivative security is decomposed in these basis instruments. The importance of these basis instruments tools cannot 
be overemphasized. A top-level understanding of their relevance may start by asking the two following questions with 
no apparent bearing on the issue at hand: 

1. Why mapping the human genome in 2000 was hailed by the medical research community as one of the major 
scientific achievements of the century in the field? 

2. Why was Mcndclcicv's (1869) description of the periodic table of elements such a foundational achievement in 
chemistry and physics? 

In answer to the first question, we have been explained that identification of the role of each gene in our DNA will help 
understand the causes of most genetic diseases, facilitating tailor made cures. Although in most cases, these benefits 
are yet to be seen, governments and private sector have and continue to invest massively in this endeavor. 

In the second case we now know for sure that because Mcndclcicv's table helped scientists understand that all phys- 
ical matter were just combinations of his basic elements, it became clear that all desired materials could be built by 
doing just such combinations. The issue shifted to economics: how to make such combinations most cost effectively? 
In cases where cost effective processes have been found, this has led to the development of chemical industries who 
have taken advantage of the opportunity and helped build the wealth of nations by satisfying our needs for all kinds 
of materials. Research still continues to expand to find more cost cfFcctivc processes to manufacture needed end products. 

In financial derivatives risk management, our identification of the basis instruments and the characteristics of each one 
of them can be likened to the identification of the gene as the unit clement in the expression of each living being de- 
scriptive feature as well as the inventory of all possiblo genes. In chemistry or physics, the analogy is the atom and the 
inventory of all atoms given by the Mcndclcicv table. Our decomposition formula is analogous to describing the genetic 
composition of each possible living being, once identified; in physics or chemistry, our decomposition formula could be 
equivalent to providing the atomic composition of each described material, whether solid, liquid or gas. 

These analogies arc intended to relate to the importance of the proposed innovation. A survey of published literature 
will help put this innovation in the proper context. 

Published literature directly relevant to Definition 1 and Theorem 1 discussed later, arc the contributions of Kenneth 
Arrow (1953, Nobel Prize 1972) and Gerard Dcbrcu (1959, Nobel Prize 1983). Arrow and Dcbrcu defined what arc 
now called Arrow-Dcbrcu primitive securities (and Arrow-Dcbrcu state prices) as the basis elements in which all other 
contingent claims could be decomposed into in a single security, single future period market (but with multiple possible 
states). Arrow Dcbrcu securities arc simply securities that pay $1 at the end of the period if a given state occurs and 
nothing if any other state occurs. If there is one Arrow-Dcbrcu security for all possible future states, then a claim 



WO 03/107137 



8 



PCT/US03/19179 



contingent on future states can be decomposed into such securities. To sec how this happens, one can observe that if f 
is a function defined on £, set of all possible states the variables S can take, we have the trivial identity: 

/(s) = £/(*) w 

where l^s=a> represents the function taking the value 1 if 5 is equal to s and the value 0 if S is equal to any other value. 

The application, which to each payoff at a future time T that is a function / of S stochastic variable known only at 
T, associates the current value of the contract is an homomorphism of group if state prices arc taken as fixed. In the 
traditional risk neutral pricing framework, this can be seen in the pricing by taking expectations and the discount factor 
J9(*o>*o) which represents the value at to of the security paying 1 with certainty at time to- We have: 

E(f(St l )B(t 0) t 1 )) = £ x l {5tl=a> B(t 0 ,t 1 )) = xPr^(fti =*)B(*h*i) (2) 

In a single period market, since every contingent claim on the value St t of the underlying S at time to will have a 
payoff that can be mathematically written at time to as a function of 5t lt i.c f(St x ), identity (1) demonstrates the 
Arrow-Dcbrcu result that any contingent claim, which can always be described by a function /, can be obtained as a 
sum of f(s) securities paying 1 at time to if St t is equal to a particular value s and zero if any other value is taken, with 
a spanning the whole range of value S may possibly take at time to- Let's further detail this concept in a simple example 

Lot's suppose we arc on Oct 27,2001. We consider the simplest case of a ticket purchased at 10AM EST(Eastcrn Stan- 
dard Time), October 27, 2001=t o representing a bet that the New York Yankees (AL) will have won the 97th Baseball 
World Scries against the Arizona Diamondbacks (NL) at the end of the day on November 4th, 2001. If this event 
occurs the ticket holder will receive $1.00 at 10AM EST on Monday Nov 5, 2001=t 7 , otherwise he will receive nothing. 
There arc only two possible outcomes/states; with probability 7r 0f i state 1 of nature, that is AL winning will be realized 
and with probability 71*0,2 state 2 of nature, NL winning will be realized. We also assume fractional dollar bets can be 
purchased. Using daily fed funds rates reported by the federal reserve, we obtain B(to,t 7 ) = 0.9994. 

So, the NY Yankee victory bet is an Arrow-Dcbrcu security in this state space and we denote it . Let's supposo there 
is market consensus on 10/27/2001 that there is a 40% (= 7r 0 ,i) chance that the NY Yankee will win and a 60% chance 
that the Arizona Diamondbacks will win (= ?ro,2). 

Since there arc only two possible states for the lottery ticket at time t 7 , there remains only one Arrow-Dcbrcu security 
to be described to have what is called a complete market, i.c. one in which all other contingent claims arc mere linear 
combinations of elements of the set of primitive securities. This remaining security is Instate 2}- 

Sincc B(to,t?) is traded in the market, i.c. deposit accounts paying the fed funds rate exist, Instate 2> is redundant 
since: 

l^State 2} = 1 — Instate 1> ( 3 ) 

Since $1 can be safely replicated on 11/05/01 with $ 0.9994 on 10/27/01, this allows a static replication of on 11/05/01 
by scaling on 10/27/01 a bet ticket for S 0.4 0.9994=80.39976 adding $0.59964 to it and depositing the sum equal to 
$0.9994 in a safe account paying the fed funds rate. Indeed if the Diamondbacks win, the bet ticket will expire worthless 
on 11/05/01 so that one would receive $ 1 from the deposit paying the fed funds rate; if however, the Yankees win, 
the $1 taken from the deposit account will be used to pay the holder of the bet ticket sold on 10/27/01. As such the 
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position built is exactly equivalent to buying a ticket paying $1 on 11/05/01 if the Diamondbacks win and zero otherwise. 

Let's suppose now that we have a derivative security paying f(St 7 ) at time tj. Since there arc only two possible states, 
/ is a binary (Boolean) function and we note fj = /(Sia*el), / 2 = f(State2). Using the Arrow decomposition, we 
postulate that holding the derivative security paying f(S i7 ) at time t x is equivalent to buying (/i - / 2 ) identical lottery 
tickets paying 1 if AL wins and putting $ 0.9994 / 2 in a secure deposit account paying the overnight fed fund rate of 
interest. Indeed, it is easy to sec that if state 1 occurs at t u one receives $ (/i — / 2 ) for having the winning tickets and 
adding to that the $ / 2 withdrawn from the secure deposit account, one ends up with $ f\. If however State 2 occurs, 
one receives nothing for the lottery ticket and ends up only with f 2 withdrawn from the secure deposit account. Thus 
the equivalence with the derivatives security paying f(St 7 ) at time £7 is evidenced. 

The problem we faced to begin with was that, since the end of the fifties and beginning of the sixties when those early 
works were being pioneered, little has been done to extend these Arrow-Dcbrcu securities as primitive securities in a 
multi-period market where there may be many securities without making questionable assumptions: 

In the early 1970's, Black, Scholcs (Nobel prize 1997) and Morton (Nobel prize 1997) made the assumption that the 
dynamic of securities followed a pattern similar to that of a stochastic process called a geometric Brownian motion with 
a constant parameter called the volatility. They thus derived what is nowadays known as the Black-Scholcs-Mcrton 
formula. What made this model appealing was that if its assumption held, the basis instrument from which to replicate 
any contingent claim on an underlying would be the underlying itself. Black, Scholcs. and Mcrton showed that a dynamic 
hedging strategy would be instantaneously risk free if it consisted of selling an option (or any other contingent claim) 
and buying the underlying in quantities equal to the derivative of the price of the option with respect to the underlying 
(also called delta) according to the formula they calculated. Under a most general extension of their work, the argument 
goes as follows: 

Let's suppose the underlying follows a diffusion process, i.e. 

dS t = v(Sut)db + a(St, t) dW t (4) 

Then if one considers a contingent claim on the security S whose price at any time t can be obtained as a function of the 
value of the underlying S at time t and t itself, i.e. /(SM)> then, if also assume / is twice diflcrcntiablc as a function 
of St and diffcrcntiablc as a function of t, then one can apply the well known It lemma which gives: 

can be rewritten, using (4) as: 

*-8S*+(5S + j'*«-'»S)* 

which means : 

= S,t + Ai) = f(S t =F) + ^(S t = F)(S - F) + (§£ + \<? (Su 0 0) Ai + 0(At) 

or : 

/(S.+A* = S,t + At) - f(S t =F)- |£(S t = F)(S -F)=(ZL + 1<t* (F, t) 0) At + 0(At) (8) 



(6) 
(7) 
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that is: 



E [(/(W = S t + At) - m = F, t) - ft (ft = F)(S - F) -(% + $««(**)&)**)'] ^ 

Urn — — T* " ~~ ® 

At— o At 

Formula (6) shows that a portfolio consisting of a long position in the derivative contract whose value is equal to f(S t ) 
at any time t coupled with a short position in f| units of the underlying is solely a function of time, thus risk-free. 
Thus it should earn the risk-free rate of return ra(i). So if one writes 

U t (St) = /(*, S t ) - j^S then <m t (S t ) = n t (S t )r d (t)dt (9) 

Since St pays a dividend at a rate named r f (t) assumed to be constant between t and t+dt, rj(t)§^dt should be 
subtracted from the variation of U t (St) coming from (6) .Thus, 

(If -«c*>(/<**>-g*) (10 > 

In the case of a Call Option struck at K and maturing at T, it is known that 

/(Sr,T) = CoXl{S T ,K,T y T) = (S T - K)+ (11) 

If it were assumed that the underlying followed a geometric Brownian motion with constant risk-free rate r and volatility 
and no drift, one would have: 

^gqlKA, Jr,t,p + l ^CaU(S t JC,t,T) j = y ^ Call(SMT) _ ^ BGaUj^M^ (u) 

(12) and (10) were obtained by Black, Scholca & Morton. Recognizing in (12) an equation reducible to the heat equation 
well known in Physics, they were able to derive the celebrated closed form formula carrying their name: 



Call(S 0> K, 0, T) = S 0 JV(d 0 + °VT) - Ka~ rT N{d 0 ) 
■with do= ^ (S "^ r -^ )T 
In obtaining this result, the following assumptions were made: 

(i) It was possible to trade in continuous time. 

(ii) Market on the underlying was perfectly frictionlcss; by frictionlcss we mean no transaction costs, no bid/ask 
spreads, no restrictions on trade (legal or otherwise) such as margin requirements or short sale restrictions, and no 
taxes. There is infinite liquidity. 

(iii) The dynamic of the underlying follows a geometric Brownian motion with constant risk-free rate of return and 
constant volatility. 

Assumption (i) has proven impossible to realize in practice because of physical limitations and the importance of trans- 
action costs; derivatives portfolio managers tend to rebalance their portfolio to adjust the sensitivity of their portfolio 
to spot changes once or twice a day. 

Assumption (ii) is just not true in practice and virtually all data scries for undcrlyings traded in open markets display 
skew and excess kurtosis incompatible with the geometric Brownian motion assumption of (iii). The failure of portfolio 
insurance in the 1987 stock market plunge, and to some extent the "normal" assumptions behind the failure of LTCM 
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in 1998 ([46]) reveal that adhering to such assumptions can have a potentially devastating effect on the stability of 
the U.S. financial system. The resulting appearance in 1987 of smile and term structure of implied volatility quotes 
in vanilla options markets and the increased liquidity of vanilla options led io the development of a model by Dupirc, 
Dcrman & Kani, Rubinstein ([25] [53]). This model is based on the assumption that vanilla options (Call/Puts) for 
all maturities and strikes arc the basis securities from which the price of all subsequent derivatives instruments, path 
dependent or not, might be inferred. This model, while having theoretical appeal, has not been validated by empirical 
tests ([24]). In fact, practitioners price path dependent options with concerns for the Vega-convexity due to the stochas- 
ticity /uncertainty of volatility in a manner that cannot be reflected in Dupirc & al, ([25])). A close analysis of the 
Dupirc-Dcrman-Kani-Rubinstcin model, if we position the problem in finite future states that correspond to the reality, 
shows that accepting their model would be equivalent to stating that a system of sn unknown would be determined by ns 
equations without any empirical or economic justification on the rule used to determine the remaining sn-ns equations. 
The most advanced models used in practice arc now going back to and extending models with Poissonian jumps or 
stochastic volatility earlier pioneered by Morton [50] and Hull and White [40]. However, Das and Sundaram [19] have 
shown that the term structure patterns of these models arc fundamentally inconsistent with those observed in the data. 
Stochastic volatility models, which fare better than jump models in data description according to Das & Sundaram, 
arc found to lead to greater incompleteness as quantified by Bcrtsimas, Kogan & Lo [6]. It now appears there arc a 
wide range of new models of markets dynamics to replace Black-Scholcs, all without clear undisputed winner for all 
circumstances. The proliferation of these models has created additional difficulties to mark non-standard derivatives 
to market, making it difficult to strictly abide by FAS 133 & 138 ([28], [30]) for accounting for derivatives positions in 
financial statements; the Financial Accounting Standards Board (FASB) introduced on June 15, 2000 [30] mandatory 
compliance with FAS 133, the latest accounting rules for derivatives to emerge since 1984. Reporting non-standard 
derivatives in balance sheets and derivatives P&L on income statements at fair market value, as is requested in FAS 133, 
can be a daunting challenge due to lack of liquidity on these products and the lack of standardization of mark-to-market 
methods. In the derivatives market the recent loss by Allfirst of $691 million in foreign-exchange losses went unnoticed 
over a couple of years due to inappropriate mark-to-markct of its derivatives positions. The existence of decomposition 
in basis instruments and the existence of a market for these basis instruments, as made possible by the present inven- 
tion, solve the problem in a definitive way. Furthermore, FAS 133 requires that derivatives be accounted as hedges and 
their profit and loss (P&L) reported in the Other Comprehensive Income (OCI) section of financial statements only 
when it can be shown that this P&L offsets that of an existing underlying in at least 80% of the risk of fluctuations 
in the value of this underlying. This creates potential confusion as to how to estimate this 80% value among the wide 
choice of possible models. Patents applications WO 02/44847 A2 or US 2002/0107774 Al arc one of the many model 
dependent approaches to address the issue. The existence of a market for basis instruments would conclusively ad- 
dress this issue. Finally, all those models, by assuming continuous-time processes arc only approximations to physically 
realizable phenomena, due to the existence of transaction costs. In fact, Bcrtsimas, Kogan & Lo's [6] quantification 
of this approximation show that even if lognormal assumptions were true, the approximation error may be of signif- 
icant magnitude for some payoff types, further reducing the effectiveness of all the models surveyed for hedging purposes. 

Moreover, in obtaining tho PDB (12) or (10), two important assumptions were made which arc not always true in 
practice: 

• It was assumed that the derivatives contract could be put in the form f(S t> t) 

• It was assumed that the asset followed a diffusion process, thus we could apply It's lemma. 



WO 03/107137 



12 



PCT/US03/19179 



la practice, for many important cases, the derivatives contracts price cannot be obtained as /(St, t) for all times t. 
Additionally, extensive mathematical finesse is required to derive an applicable PDE for each new derivatives contract 
under the diffusi on process assumption. Such extensive mathematical methods can not be not readily automated. The 
approach to obtain the PDE for the Average rate or the Lookback option may exemplify this difficulty. Sec for example 
the PDEs compiled in Rogers & Talay [52]. When no diffusion process assumption is made it is substantially more diffi- 
cult to obtain a tractable PDE. One must then resort to trees or Monte-Carlo computations for pricing purposes. More 
significantly, it is difficult using these methods to understand how to incorporate micro-structural issues in the pricing 
of derivative contracts. Our invention docs enable easy incorporation of micro-structural issues when pricing derivatives. 

The current approach for managing the risk of a portfolio of derivatives is by hedging the Greeks of the portfolio. The 
Greeks arc the sensitivities of the portfolio to different market variables. Delta is the sensitivity with respect to variations 
in the asset price. Gamma is the sensitivity of the Delta to variations in the asset price. Vega is the sensitivity of the 
portfolio to variations in its Black-Scholcs-Mcrton implied volatility. Thcta is the sensitivity to time. Higher order sen- 
sitivity parameters also exist, with non-standard names. In general, hedging strategics arc aimed at matching first and 
second order derivatives of portfolios to the various parameters. In doing so, what is implicitly assumed is that the sensi- 
tivity of portfolios of derivatives to their parameters is polynomial and of degree 2 or 3 or 4 for the most accurate hedges. 

For example, one would take a portfolio n(S f) tf) and write 
II(5,t,a,r) »n(S 0 ,tn,cr 0 ,r 0 )+ 

< ^ flnfgo^o.ao^o) ^ aruso.to^.ro) ^ ancso.to.cro.Tm) ^ en(So,to,<r 0 ,ro) ^ > (5-. J g 0> t-t 0> <j — <r 0> r-r 0 ) > 

+* (S - So, t - to t a - cr 0 , r - r Q )H (S 0i to, a 0 , r 0 ) (S - So, t - to, a - cr 0 , r - r 0 ) 4- ... 

where H is the Hessian of n with respect to the variables S ti <r, r. 
and then one would consider that a perfect hedge for n(S, i, s) would be a portfolio matching the value at the initial 
point as well as the first and second order sensitivities. In reality however, the hodgc is dependent on the model used to 
infer the dependence on the specified parameters. Additionally, the polynomial approximation derived from the Taylor 
expansion is of uncertain accuracy and can lead to serious problems, especially for highly non- linear pay-off options such 
as barrier options. Preferred embodiments of our innovation address this problem as well. 

The idea of using basis elements or spanning for the pricing of derivatives continues to this date to be an active field of 
research in Financial Economics. For examples of recently published approaches sec Madan, Carr, Gcman, Yor, Bakshi 
[48], [13]. However, by making assumptions of continuous time hedging, brownian filtrations or semi-martingale prop- 
erties of the underlying and absence of micro structural effects, these approaches not only limit the possible accuracy 
obtained for pricing purposes, but they create a model risk when one would attempt to use those results for hedging 
purposes using the selected basis (Characteristic function or Hcrmitc polynomials as in [47], [48] for instance) 

Our invention differs substantially from the approaches mentioned in the prior art papers by the fact that wc focus 
on selecting exactly replicating basis instruments that arc as close as possible to what exists on the market, i.e. spot, 
forwards/futures and vanilla options. Additionally, our implementation is designed for computational tractability. In 
the prior art papers, the authors make simplifying assumptions for mathematical tractability that arc not always eco- 
nomically justifiable. A recurrent feature of such decomposition approaches is to seek Hilbcrtian bases and obtain 
pricing by taking the projection on a few elements of the basis without clear justification of the economic rationale of 
the approximations. Furthermore in those approaches, one has first to non- trivially infer or imply the price of basis 
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instruments from traded instruments. 

In their simplest form in a two state, two periods economy, a preferred embodiment of the basis instruments we introduce 
could be cquivalcntly defined in addition to the Arrow-Dcbrcu primitive security as: 

An agreement at time t 0 to buy at time t x a lottery ticket at the price $10,000.tti,i if state 1 occurs and $10, 000.7Ti, 2 */ 
state 2 occurs at time t\. 

When extending to n arbitrary periods, a basis instrument would be: an agreement at time to to buy at time tk-i 
(1 < k < n) a lottery ticket at the price $10, OOOttA: - 1, f k -i, where f k -i is a k-1 series of 1 and 2 indicating the states 
. of winning or losing lottery tickets bought between time to and time tk-i . 

How far into the history of these tickets wins or losses the seller will take into consideration will depend on his views of 
how past events affect the likelihood of future events. 

Our example of the 97th baseball World Scries can be used to further illustrate the concept. To win the World Scries, 
a team needs to win four of possibly seven games. For the 97th scries, the games were scheduled as: 



Game 


Date 


Location 


Game 1 


Sat, Oct 27 


Arizona 


Game 2 


Sun, Oct 28 


Arizona 


Game 3 


Tuc, Oct 30 


NY Yankees 


Game 4 


Wed, Oct 31 


NY Yankees 


Game 5 


Thu, Nov 1 


NY Yankees 


Game 6 


Sat, Nov 3 


Arizona 


Game 7 


Sun, Nov 4 


Arizona 



Table 1: World Scries 



We call Xt the stochastic process taking the value 1 the day after the end of each game won by the NY Yankees (AL) 
and the value 0 otherwise. As such X t is for * being a day after a game a bet that the Yankee will win the game of the 
day before. State 1 that AL wins the World Scries is thus such that: 

l<Stotei> = ^Ej^-xy,^} — /(-Xti,-Xta»-».-X* T ) ( 15 ) 

which shows that the Yankee World Scries victory bet is now a path dependent derivative contract on the realizations 
of the stochastic process X t * Let's suppose that to replicate this contract the only instruments we have arc bets on each 
individual game and that those bets arc taken at time to =Oct 27, 2001 at 10AM EST (Eastern Standard Time) but 
paid for one period before the game settlement date. It is natural to expect that as games arc played, estimates of the 
probabilities of victory in the subsequent games change. 
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1.1 Notations 

For i = 0..6 wc note: tt* = 7r j>i+ i(Xt i+x = l|^t,-, X to ) so that 

Kit+liXti^Xu, ..^to) = TTi.i+l = Xti^lTi + (1 - - TTi) 

i.C: TT^i+iCX^J = (1 " TTi) + (27Ti - l)Xt, +1 

For any derivative contract on the stochastic process X t paying at time £7 f(X tl ..Xt 7 ) , wc note n{ the price at the 
time U of the contract paying £7 f(X tl ..X t7 ) at time £7. 

Wc suppose each team has a 60% vs. 40% probability of winning its first match in its home stadium and that this 
probability is increased by 10% after each game won if the next game is at home and it decreases by 10% after each 
game lost at home if the next game is also at home. Wc also assume that once a party wins 4 games, its probability of 
winning the remaining games if any is zero. These translate algebraically into: 

7T 0 = 0.4 
7T 2 = 0.6 

For i — \ and 3: 

TTi = (0.2X ti -0.1)+7Ti_l 

For i = 4 and 6: 

tt< = {{Q2X U - 0.1) + x l {i - 3 <z;}»i^ J <4> + l 4*-3>E}oi** ) 
Now wc introduce our 2-statcs multi-periods basis instruments definition. 

1.2 Definition 

In a single process (X t ), two states (1 and 0), multi-periods market {£0 < £1 < < ^} with reference currency A?, 
so that X? is equivalent to 1, a basis instrument is defined as: 

An Agreement between two parties, 0 and £1 contracted at time £0 and stipulating that cither: 

• At time U-i, 0 shall pay Q the amount N(X tl , ® *£i(X tl i ...,X tM ) and at time t it CI shall pay 0 the 
quantity : N(X tl , X^s) x X^ 

or 

• At time ti-u 0 stall pay Q, the amount Tf{Xt t *... % Xu^ x )& *i-i(X ti , ~,X Umml ) and at time £», SI shall pay 0 the 
quantity : 7?(X tli ....^..J x X° 

When Q pays 0 the quantity NiX^ , ... , X tM ) , it is an extension on what is commonly know as a forward rate agreement 
and will be referred to here as a zcroth order basis instrument in the reference currency 
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When Q pays 0 the quantity N(X tl , ..^Xt^), it is an extension on what is commonly know as an Arrow-Dcbrcu 
security and will be referred to here as a zcroth order basis instrument in the currency X}. 

1.3 Proposition 

In a single process (X t ), two states (1 and 0), multi-periods market {to < *i < *n-i < wc navc tnc following 
relationships: 

And for 1 < i < n, 

4-1 = (4(*ti> 1) - *£(*tn ».-Xi.-i.O)) 8> *& l + n£(X tl , ..,Xt^, 0) 0 Trfjr 1 

If wc apply this to our World Scries example, to calculate the price at tO of an AL world scries victory bet, wc have: 

f(X tli Xt 2i ...,Xt r ) — Iffij^x*^} 

and using our formula (16) with data (15) wc obtain a probability of 48.3% and therefore if there is a market where bets 
on the AL victory exists, a risk free arbitrage would be to buy a ticket betting an AL victory at 40% of the face value sell- 
ing basis instruments replicating that victory bet at 48.3% of the face value, realizing a net profit of 8.3% of the face value. 

The example above presents the concept of basis instruments and how this innovation is used for pricing and hedging 
any derivative security in the case of a single process, two states (1 and 0), multi-periods market in a way that is an 
innovation on the Arrow-Dcbrcu primitive securities. To infer conditional probabilities of possible future occurrences 
from basis instruments prices, our new method is to use the relationship 

Prob(X t , = l|*t<>> - >^"t«-i) = go 1 

In the more general case of a multi-processes, multi-states, multi-periods market, in another embodiment of the inven- 
tion, wc will present below the more general innovation that applies. 

While in a two state process, vanilla options and Arrow-Dcbrcu primitive securities coincide, when the number of states 
increases, the two types of contracts become different, yet the two, with bonds and forwards added to options, remain 
equivalent through a bijective transformation. Wc choose to present our basis instruments below as extensions on 
vanilla European options forwards and cash instruments instead of the academic Arrow-Dcbrcu securities as this basis, 
is more economically sensible and computationally allows a faster convergence to the desired payoff for most derivatives 
instruments using the lower order basis instruments and making the decomposition around the forward value around 
which probabilities of occurrences arc clustered. This is also of critical importance to relate to the shortcomings of 
the Black Scholcs Morton derivatives hedging approach. Formula (14) shows that at any trading time, a derivatives 
position would be hedged through cash and underlying or forwards in amounts equal to the sensitivity of the option 
to its underlying (or delta), but unlike Morton, there arc complementary terms representing a spanning of calls and 
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put options in proportions equal to the second derivative (or gamma) of the financial derivative with respect to the 
underlying. Formulas resembling (14) have been known before in a two period case, (sec for example [16], appendix), 
however these formulas have been presented with two important limitations solved in our derivation: 

• It is assumed that the derivatives terminal payoff as a function of the underlying asset is twice diffcrcntiablc 
everywhere, substantially reducing the number of real life derivatives to which the formula can be applied. It is 
also assumed that the asset values state space is continuous. 

• The formula is derived solely in a single asset, two period setting and thus merely presents variations around 
Arrow-Dcbrcu. 

2 General notations and definitions 

We suppose we arc in an m+1 assets economy, with n+1 trading periods to < < U ~ < *n chosen as seen appro- 
priate. 5 C taken by convention to be So is the base currency asset and thus merely represents the unit of cash in the 
reference currency and S 1 , ...,5"* arc risky underlying?* whose values in units of So vary across time; so, the default unit 
of any given number or quantity is So and when there is no ambiguity, Sq is simply one. 

JSb.n represents tho m+1 undcrlyings economy between to and ^. e 0 ,n represents an exchange system of contracts whose 
values is dependent on the realization of the m+1 imdcrlyings of economy Eo t n between to and tn. 

refers to the value of the underlying S i at time tj. S) refers to the realized value of the underlying S) at time tj as 
opposed to the abstract reference to the parameter value symbolized by 5} . 

Flj is an arbitrarily chosen value of the underlying S* at time tj. In general for the developments to be made here, this 
value wiD tend to be the forward value for the maturity tj at time or the value of SJ.j. 

8 = (Sj) x<i<m S = (§]) 1<i<m F = {Fj) i<i<™ S 3 = (§}) 1<i<m Fj = (Fj) !<<< m 

A"+ is the maximum of the real value X and zero, or Max(X t 0). If / is a function denned on [a, 6], then the real 
b 

quantities noted as / f(x)dx and §£(2*) arc subordinated to the definition of a finite set 7 X = *o = o. < xi < < x H = b 

a 

so that ff(x)dx is by definition equal to: £ /(x^ x (x i+ i-Xi) and |^(x £ ) is by definition equal to ^^jlff** ; s mav 
be equal to the infinity and will be truncated to a finite number only upon proof that the remaining terms arc negligible 
within the approximation bounds sought. In general for the developments to be made here, Ix will tend to represent 
the actual range of values an underlying can take or any bijective transformation of this range. Since increments of an 
underlying valuo for all types of undcrlyings to be considered arc discrete (that is there is a minimum increment value 
called the basis point), the cardinal of I* is always finite or dcnumcrablc ; 

/ /(x)dx may also be formally noted as / f(x)dx\ Tm . H(x») may also be formally noted as %(xi)\ ItB For any x G I x , 

tt o 

there exists a j < n 

Multiple derivations or integration signs for functions of several variables simply mean that the operation is iterated on 
the specified variable in the integration or derivation symbols. / .. / /(x 1 , .., x k )dx 1 dx k may also be formally noted as 
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/ . . . J /(x\ .., x^dx 1 . . . dx*\ Iml , .. tJ> . 

inay 3130 DC formally noted as ^r(^,..,^)|/ s i , 
Thc indication / b1 , ..,/ sP on the integral may be omitted when the subordinated intervals partition is evident in the 
context or irrelevant. If x clement of I x is equal to x,-, then x + is by definition Xj+i, X++ is by definition Xj+2 and x_ 
is by definition Xj_i. If j = s, x + is by definition x^-, x++ is by definition x 3 -; if j = 0 x_ is by definition Xj. 

If p is a real, sign{p) = + if p > 0, sx^n(p) = - if p < 0, si^(p) = (void), if p = 0. The convention x«^( P ) is used in 
the result presented below. = between two quantities identifies an equality between the two quantities which is true by 
definition. For any set A included in a larger set X % for any x e X , we define the function 1^ as 1a 0*0 — 1 if x € A , 
(x) = 0 if not. If A is defined by a Boolean condition, A in I a may be replaced by that condition. 

The key to the solutions described in the present invention as will be further detailed below is the introduction of 
functions of real numbers rather than just real numbers in expressing payment terms in premiums of contracts in a way 
that has practical m eaning . 

3 Basis Securities definition and notations 

In order to make the d efini tion more accessible wo stress characteristic features of the Basis instruments as used in a 
preferred embodiment of the present invention and the rationale for those features; 

(1) A Basis instrument Contract (BIC) is a contract and involves and identifies two parties named buyer /? and seller Q 

(2) The definition of each contract comprises three dates: 

• The contract agreement date to, which is the date at which the binding rights and obligations on both sides of 
the contract arc agreed upon. 

• The premium payment date U , with U > to , which is the date at which the party identified in the contract as 
the buyer /? complies with its part of the agreement by paying the seller Q an amount in units of basis currency 
known as the premium of the contract. 

• The contract expiry date also known as maturity tj, with tj > U , which is the date at which the party identified 
in the contract as the seller Q complies with its part of the agreement by paying the buyer P an amount in 
units of basis currency defined as one of the generic form further detailed below. 

This facilitates the formation of BICs between trading parties and can be implemented in a system or computer program 

product. 

Definition 1: 

A basis instrument Bf 0tUttj (ft Q, N, (ij , . . • , i k ), (S x , • • • , S k )j (K u ■ • ■ , ^fc)» (h » * * ' » J'p)) fe a sccurity.contract entered into 
at time to between a buyer /? and a seller Q stipulating that: 

• At time i*, /? shall pay CI a premium payment amount such as , 

Bt M (ft ft, iV, (n, . . . , i k ), (6 1 , . . . , 6 k ), (K u • • • , K k ), (j u • • • , 3p)) (ft; - ■ ■ . *, 

units of underlying currency S c . 
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• At tdmc i», Q shall pay P a payout payment amount such as, 

units of underlying currency S c And where <p , called the scaling density function is an application such that: 

Bio.ufy W, (l,N,(ii,- ',ik),(Su - ,Sk),(Ki, - ,K k ),Uu- -, 3 P )) (ft; • • • . ft) m 

N(S tl ; - -, ft) ® Bf BtU>tj ((»„•■•, t k ), (*„ • • • , ft), fri. ■- , K k ), Uu • • • . Jp)) (ft: — . ft)S£ = 

0 

x**W; • * " • **)■ C*. • ' • . «b). ■ ■ ' , Kk), ft. • • • , Jp)) (ft; • ■ • , ftJSg 

Where N 0 * as the inventory of the of the counterparty selling 

BtG % u ,t, (A ^» N> ft , • • • , i k ) , ft , • • • , 6 k ), (Id , • • • , K k ) , ft , • • • , j p )) prior to the said transaction. 

^Sia a (&> ■ • ■ i (*n- ■ • . *)» • ■ ■ i **). (iw • • .ip)) (ft; ••■,*.) 
= i£ fi^,*,*, ((n, • • • • • ■ . <W> • • - .**). 0"i» ■ • ■ , i P )) (ft; • • • . ft<) 

The set of Basis Instrument Contracts is called a BIC -Basis. 
When the payout payment amount is in the format 

n (aft,. • • ,flfc, • .. . 8%, ■ ■ ■ ,sfr) (ft (fl* - Jfi)) + • • • (ft (flfc - #*)) + s£ . . . s£sg 

, then we term it an Extended Option Format Basis Instrument Contract Payout or EOFBICP. It naturally imply the 
premium payout format 

(P> n.*. (*i> • • • . *0. ■ • • . 4). • ■ ■ . **)> Oil • • • .*)) (ft; • • • . s u)si 

= N (S U • • • , ft,) (ft,. • • , u), ft, - . - , fa), (JKi, - - ■ , Jf*), ft, - ■ • ,i P )) (ft; — , Su)St 

The EOFBICP is used here to facilitate teaching of the present invention. Other payout formats can be used to create 
a BIC-Basis. Such format may include Extended Arrow-Dcbrcu Format Basis Instrument Contract Payout or EADF- 
BICP, the Extended Fourier Transform Format Basis Instrument Contract Payout or EFTFBICP, the Extended 
Hcrmitc Polynomial Format Basis Instrument Contract Payout or EHPFBICP. All those format arc equivalent and 
from a complete set of one kind, a complete set of the other kind can be obtained through a linear transformation. In a 
preferred embodiment of the present invention a few examples of the linear transformation arc disclosed in a matricidal 
form. 

The scaling density function translates how the offer and demand for each basis instnuncnt in turn affects the level of 
prices. This phenomenon is also known in securities markets as slippage. For 0 < p < & + 1, if p is a priori a function of 
(S*i5 ••"!&«)• 

& (ftiJ * ' * » fti) k called the notional of the contract. 

W (ftn • • ■ » ft<) & ((*!»' • ' • ik), ft, - • ,5k), (Ki, . - • , K k ), ft, ■ • • » 3 P )) (ft; ' ' • » ^ the premium or price 

of the contract. 

fc is the order of the basis instrument. When k = 0, the basis instrument is simply noted as B^ tUftj (N) and when 
in addition, 

V> a* • ■ <*). (*i , — . 0i, • * • .*))) (JV (Si ; - - , S u )) = 

JV(ft; ■ • - B^ >ti>ti ((n, . . . t4fc ), ft, . . . , 6k),(K u -*-,Kk),tii, --J P ))(S 1 ; - ,St i ) 
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, the basis instrument is simply noted as B^ tUtti . 

When the basis instrument price at U does not depend on (S\\ • • • , St<) as is the case for zcroth order basis instruments, 
the notation of the price is the same as that of the basis instrument contract. 

In current markets, for most market makers, <p is a stepwise increasing function, —<p(—N) — <p(N) > 0 and is called the 
bid/offer spread. For practical implementation purposes, in a preferred embodiment of the present invention, may be 
defined and input in an implicit or explicit manner. Examples of an implied input or definition arc provided further below. 

Recent patents have addressed the problem of financial market liquidity in the context of an electronic order-matching 
systems . Patents relating to derivatives, such as U.S. Pat. No. 4,903,201, disclose an electronic adaptation of cur- 
rent open-outcry or order matching exchanges for the trading of futures is disclosed. Another patent, U.S. Pat. No. 
5,806,048, relates to the creation of open-end mutual fund derivative securities to provide enhanced liquidity and im- 
proved availability of information affecting pricing. This patent, however, docs not contemplate an electronic derivatives 
exchange which requires the traditional hedging or replicating portfolio approach to synthesizing the financial deriva- 
tives. Similarly, U.S. Pat. No. 5,794,207 proposes an electronic means of matching buyers' bids and sellers' offers, 
but outside the scope of an electronic derivatives exchange which requires the traditional hedging or replicating port- 
folio approach and without explaining the nature of the economic price equilibria achieved through such a market process. 

U.S. Pat. No. 5,845,266 and U.S. Pat. No. 6,098,051 implement order-matching and limit order book algorithms, which 
can be and arc effectively employed in traditional "brick and mortar" exchanges. Their electronic implementation, 
however, primarily serves to save on transportation and telecommunication charges. No fundamental change is contem- 
plated to market structure for which an electronic network may be essential. Second, the disclosed techniques appear 
to enhance liquidity at the expense of placing large informational burdens on the traders (by soliciting preferences, for 
example, over an entire price-quantity demand curve) and by introducing uncertainty as to the exact price at which a 
trade has been transacted or is "filled." On the contrary, a preferred embodiment of the present invention is to provide 
through the scaling density functions a functional format that reduces the size of the data to be supplied to the system 
and in effect eliminates the need for a trader to continually update their bids and/or offer responsive to demand. In 
fact, once the scaling density function has been provided for all basis instruments contracts, there is no need for further 
trader intervention and the trader or market maker knows very precisely for each quantity transacted at which price the 
transaction concluded. 

What is different between this type of contract and other types of derivatives contracts is the functional nature of the 
premium as expressed in the definition of the contract and the introduction of the date t* as separate and posterior 
to to. While it is true that in most OTC contracts there is a contract date and a settlement date different from the 
contract date, the difference - with aasuming as is done for modeling purposes that the contract date and the premium 
payment date arc the same, is minimal and it is straightforward to make the adjustment. So, each of the dates t» and 
tj may actually be coupled settlement dates a few days later. 

The premium is expressed as a function of the values of the fluctuating variables between time to and time U . Enabling 
this extension is what makes possible the most critical sections of this invention. Moreover, the payment at time tj 
is expressed as units of a specific basis set of functions of the fluctuating variables at time tj with scaling units being 
functions of the fluctuating variables between time to and time t». The choice of the specific functions is made so that 
the selected basis set is a mathematical basis of the vectorial space of functions of the fluctuating variables at time tj 
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with the scalar units belonging to the set of values (i.e. ) taken by functions of the fluctuating variables between time 
to and time Since in a vectorial space all bases arc equivalent, once one has been selected, any other basis would 
be equivalent to that one. As a result the claimed scope of the invention is not intended to be limited to the specific 
choice made herein for illustrative purposes, but rather to encompass its possible equivalent. In the specific case of 
the illustrative definition above, the basis is selected to be as "close" as possible to known instruments actually traded 
in the markets, that is bonds, forwards/futures and European vanilla calls or puts. The notion of "close" is meant 
in particular in the sense that if the scalar set becomes the field of real numbers, as happens at U , then wc would 
recover classical bonds, forwards /futures and European vanilla calls or puts. When tj - U = 1 day, wc would have the 
known overnight options, and overnight rcpos commonly traded in OTC markets. In doing this, wc express preference 
for hedging purposes, over alternative existing approaches where the basis is selected to be hilbcrtian with respect to a 
defined scalar product, and merely has computational appeal for pricing purposes. 

The functional nature of the premium and the scaling units is essential for the derivation of the decomposition formula 
from which the basis instruments derive their usefulness as hedging instruments. 

V 

3.1 Methods of supplying BICs 

As seen above there arc various ways of supplying BICs depending on the structure of data available or computational 
consideration made, one might prefer one over the other. Wc show hero an example of how to move between EOFBICP 
format, EADBICP format or EFTFBICP format. 

3.1.1 Correspondence between EOFBICP format and EADBICP format 

Wo have the following result, establishing the correspondence between basis instruments prices and unit underlying- 
conditional probabilities that extend and generalize the Brccdcn-Litzcnbcrgcr formula ([10])in a multi-periods, multir- 
undcrlyings setting: 



Pr oh (Si = Ki/Si-n . . . ; S 0 ) = 

TT ( K? —-if \ (16) 
^±—^6%, . . . d* K „N &Bl itM)ti ((1, . . . , m), (61 . . . , - « («« «"> 



This formula shows that with the price of basis instruments non-ncccssarily linearly dependent on the notional, the 
conditional probabilities arc also dependent on the on the notional amount. 



In a discrete environment wc further obtain matricial relationships. Wc show here how to transform conditional Arrow- 
Dcbrcu state-prices into Basis instruments prices & recovering Arrow-Dcbrcu state-prices from Basis instruments prices. 
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Wc define the conditionnal basis instruments prices vectorc as B ti _ lfti and the conditionnal state-price vector as Pa u _ uti : 



For * > /i, B? <Tljt4 (l, A:) = 0 and for fc > h, Fofc k) = 0; 

Fbr ft < i, Bg J 1(t| (1, ft) = 0 and for & < fc, Pog 1-1 A (1, ft) = 0; 

Proposition 



Wc have the matricial products: 



with: 



» 



T u l uU (U) - x^ lA (I + l)-«^ lA (I) ; 



+ 



forain + l<fc<0; Kfc</i-l 

2CtA(*.i) = ^^(ft + 1) - x ti _ 1>ti (&) ~ ~ = 1) 

3 + I" 



(fc + 1) - X tt _ ltti (fc) ^ti-Lt^fc) ~ - 1) ' 

1 3 

**l-l,tl(0)- (-1) ^-l,t 4 (l) ~ (0) 

, ^ 

(0) - **<-iA ("I) x U-uti (1) - (0) ' 

* > « (1) _ Xti _ 1>u (0) + a*.,* (2) - *,_„*,(!) 
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, * 

5? SJ . 

1 

+ ^-,,t 4 (A-l)-^-.,t < (h-2) : 

with = 1 if n = m and = 0 f or any other value of n. 
Multi-dimensional Case 

Wc indicc T and its inverse to the basis instruments considered and apply it m times to obtain the basis instruments or 
the states prices of the desired dimension. The components of B and Pa in the representation above arc m-1 dimensional 
matrices. 

Writing the compression algorithm: 

To obtain Pi(n — 1), a simple transform may help by decomposing in the basis instruments and treating the case of 
moments or baskets separately. 

After that, based on the assumption that the Pi(i) may be analytic in the dependent variables, make polynomial 
interpolation by arcs and selecting the next points based on a condition being met (for examples on the derivatives at 
the previous point). 

3.1.2 Correspondence between EOFBICP format and EADBICP format 
Let's set (ft, F t p) the probability space. 

ft™ = {So, Si, S Tt -i} 

Lot's call X the random variable with values in P. The characteristic function of X is defined as : 

= E(e«*) 



Lot's obtain all the state probabilities form the determination of To do that, wc simple assign to z, itcrativcly the 
values 0, 2tt, 2(n — 1)tt an wc obtain: 

$x(0) = l = E(l)=po+Pi + ...+p„-i 
^(27r) = E(e 2iwX ) = p 0 ^ So + pae 2 ^ + ... + Pn-i^ Sn " x 

$x(2(n - 1)tt) = ECe 2 ^- 1 ^) = p 0 e 2i < n -» s ° + p^ 2 ^" 1 ^ + ... + p n -ic 2i < n - 1)s »-* 
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This is 



( 1 \ ( 

E(e 2i ^) 



1 

p 2iwSo 



\ ( PO \ 



Pi 



V (e^r- 1 ... ( e 2^-i)«-i y \^ p^ ; 

Where wc can note: 

EX = MP (17) 

And wc can recognize in M a Vandcrmondc Matrix. For a matrix M, wc note My its cocfficcnt of line u and row v, so 
wc have here : M uv = e 2iir(«-i)S tf _i 



As wc note M" 1 the inverse of M, wc want to obtain M~ x so: 

P = M~ 1 EX 

A simple method of solution of (17) is closely related to Lagrange's polynomial interpolation formula. 



Let P v (x), v G [0, n — 1], be the polynomial of degree n — 1 defined by: 



e 2iwS v — e 2iirS 11 



Wc can define with the coefficients of the interpolation Amk, fcG [l,n],m6 [l,n] a square matrix A of dimension n. 

Wc have P„(e 2iwS «) = <f uv = E&=i A v+1>fce 2i ^«( fc - 1 ) = ELi 4i,Jfe%+i for u and v in [0, n - 1]. 
Since 6 UV = wc can replace u 4- 1 and u + 1 by u and v in [1, n]: 

n 

This says exactly that A = M _1 . Therefore, the solution of (17) is just P = AEX so wc have: 

p^^^E^ 2 ^- 1 ^) 

Application to numerical computation 



A simple algorithm to solve the problem is given in [31]: 

for k=0:n-2 

lor u=n-l:k+l 

p(u)=p(u)-exp(2i\pi S[u-l]).p(u-l) 
end 

end 

for k=«n-2:0 

for u-k+l:n-l 

p(u)=p(u)/( exp(2i\pi S [u-l] )-exp(S [u-k-2]) ) 
end 
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for u=»k:n-2 

p(u)=p(u)-p(u+l) 
end 

end 

This algorithms requires 5n 2 /2 flops. An implementation of the algorithm that yields P can be found in Mathematical 
recipes in C, Second Edition [55] : 

# include "nrutil.h 0 

void vander (double HD » double p[J, double EXD , int n) 

Solves the Vandermonde linear system. Input consists of 

the vectors M[l..n] and EX[l..n]; the vector p[i. .nj is output. 

i 

int i,j,k; 
double b,s,t,xx; 
double *c; 

c=dvector (1 ,n) ; 

if (n = 1) p[l>EX[l]; 

else { 



for (i=i;i<==n;i++) c[i]«=0.0; Initialize array. 



c[n] ** -M[l] ; Coefficients of the master polynomial 



for (i«2;i<*^;i++) { are found by recursion, 
xx « -xCi] ; 

for (j=(n+l-i);j<=(n-l);j++) c[j] += xx*c[j+l] ; 
c Cn] += xx; 



> 



for (i«=l;i<=n;i++) { Each subf actor in turn 

xx«M[i3 ; 
t=b=1.0; 
s=EX[n] ; 

for (k==n;k>«2;k — ) { is synthetically divided, 



b=c [k] +xx*b ; 



s += EX[k-l]*b; matrix-multiplied by the right-hand side, 



t=xx*t+b; 

> 

p[i]=s/t; and supplied with a denominator. 



> 

> 

free_dvector(c,l,n) ; 



> 



Int n; 



dimension of the matrix 



double* M,EX,p; 
M^vectorCltn) ; 
EX=dvector(l,n) ; 



coefficient of the Vandermonde Matrix 



input of Laplace values 
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for (int u»l;u<«n;u++H 
M[u]«exp(2i\pi S[u-1]); 
EX[u3=Value0f E[exp(2i\pi <u-I)X] ; 
p[u]«0.0; 

> 

vander (M , EX , p . n) ; computation of the probabilities . 



Regular space distribution of states 



If wc suppose that the S k arc regular distributed, for example: S k y = So + fcA. We have for the coefficient of our 
Vandcrmondc matrix: 

M uv = e 2i«"(u-l)S»-i = e 2iir(u-l)So e 2iir(ii- 1) A 

If wc define D = diag(l, e 2 ^ 50 , c 2 ^- 1 * 50 ) so that D uv = 6 uv p^~ 1)Sq . And let's have the square matrix W defined 
by W uv = with o; = e 2iTrA . Wc have: 

/ 1 ... 1 ... 1 N 



W = 



^ 1 ... W (n-1)( V ~D ... ^(ri-lXTi-l) ^ 

The u-th row and v-th column term of the product of D and W is given by: 

(DWU = Y^DukWkv = e M-i)Sa w (-i)(-D = M uv 

So wc have M — DW. If wc take further the unit step A as A = £ so that u is a n-root of unity let's note W its 
conjugate. Let's compute the term (WW) UV of row u and cohimn v of the product of these two matrix : 

n 



fc=l 



= z 



.<u-l)(w-lW*-l)(»-l) 



fc=i 



k-i 

n if « = v 

1 7l"T-T = 0 if t^t, as (««-«)" = (<r~" = 1 
= <5 ut ,n 



-I 



Aoid wc obtain: 

WW = nLd and W 1 = 

n 
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So: 

M- 1 = ^WD" 1 
n 

with D~\ = diag(l, e - 2<wSo , e -2M«-DSb). 

That gives a way to compute the state-prices p 0 > *.pn-i in only n 2 operation. But we can sec that the matrix M~ l that 
is obtained is the matrix of the Fourier Transform of order n that means that using a Fast Fourier Transform 
Algorithm, we can compute the product: 

in only 0(nlog(n)) operations. 

If we have a distribution of Sk non-rcgurlarly spaced where we write A k = S k — So but that verify the following property: 

V&€ [0,11-1], A fc GQ 

Let's set D k G Q and N k G N such as A fc = Let's set N = lcmjb e [ 0 ,„_i}(JVjb) where 1cm is the least common multiple. 

Thon we can introduce a regular division (Sk)ke[o,N-i) that fits all the states Sk and we can apply the previous result 
as we know how to compute the missing values of the characteristic function at any point. 



3.2 Incorporation of offer and demand sensitivities in the pricing, risk management and 
trading of derivatives contracts 

A preferred embodiment of the present invention solves one issue that must be addressed in the determination of basis 
instruments price. The problem is that basis instruments prices fluctuations in response to offer and demand sensitivities 
must be such that no arbitrage opportunity is created. This is tested by making sure that state prices arc positive and 
the sum of implied probabilities of all states is one. 

Although in the definition above the scaling density function is defined in an explicit form, in many practical applica- 
tions, to take into account the non arbitrage requirement the scaling density function will be formulated in an implied 
manner. It is an embodiment of the present invention to achieve this by introducing what wc define below as weighting 
functions 

Implied Definition: Weighting Functions Definitions 

Let's define a discrete probability space B = (O, B, P) with Q = 1, ,n; we define a weighting function defined on £1 x 
as W(i>rii) such that: 

Pi = 

T* 

with W{i y ni) > 0 and £ W(i,ni) > 0 to prevent negative or infinite probabilities that would create arbitrage oppor- 
tunitics. 
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Examples of Weighting Functions 

To introduce this example wc first establish the following proposition: 
Pr oposition : 

We define the target probabilities (Pi*) 1<i<n andthe base probabilities 

of themarrkctmffi)^^ for the state space under consideration 

with associated notional numbers n? 1 such that for cachi, p™ = -tH — : 

2-t 



There exists a sequence of numbers in™ , such that for eachi.pl* = 



^=1 



Tl 

and for any given value of W T » the W T ara uni Q ue ty determined for each % 
Proof : 



„ ( E ( n \ 

Si WknT — / 1 ^( ; PKe assume for instance : I 53 u;™ ) = n, 



Til' 



thus : = ^v^r" 1 ( ; »r = 



Thus for example, in a market where a first set of state priccsp^ for any first unit of Arrow Dcbrcu security is postulated, 
a weighting function may be defined as 



to reflect the fact that after state prices react ot state price responsive to offer and demand in a manner similar to the 
distribution of payouts in a parimutucl game. 

In a preferred embodiment of the present invention, weighting functions may be used to simply define scaling density 
functions of basis instruments by first transforming the input weighting functions into conditional state-prices vector 
and then, using after multiplication by a matrix T, the vector of basis instruments prices responsive to inventory, offer 
and demand for each basis instrument. 



3.3 Incorporation of credit risk in the pricing, risk management and trading of deriva- 
tives contracts 

Another embodiment of the present invention is its ability to include most accurately credit risk sensitivities in the 
pricing risk management and trading of financial derivatives. 

Credit risk in this setting is the risk that a Counterparty may not honor their financial commitments in full on due date 
and at due time. 
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A counterparty to a derivatives transaction typically assumes the risk that its counterparty will go bankrupt or not be 
able to meet its obligations at agreed upon times during the life of the derivatives contract. 

Margin requirements, credit monitoring, and other contractual devices, which may be costly, arc customarily employed 
to manage derivatives and insurance counterparty credit risk. In contrast to US Pat. No. 6,317,727 Bl, US and US Pat. 
No. 6,321,212 Bl, in a preferred embodiment of the present invention, for pricing and hedging purposes, credit risk in 
derivatives transactions is recognized as a multiplying underlying whose value in the default free reference currency and 
counterparty may fluctuate between 0 and 1, this has the benefit of not only facilitating addressing the concern of other 
methodologies, but better still, enable derivatives trading on said credit risks. 

More specifically, to take this notion into account, for each contract between two parties O and p , wc define the stochas- 
tic variable Sf ,n e [0, 1] as the percentage of a unit notional liability of N to 0 that counterparty Q meets at time t. 

In the context presented here, Sf ,n can simply be viewed as a new underlying, so that a commitment of /3 is accurately 
understood as risk free units of underlying Sf ,n . So, by multiplying the payoff f of any derivative contract by and by 
denning its scaling density function, the credit risk of the transaction is automatically included. 

So a market will fully price credit risk for all types of payoffs if basis instruments for all stochastic variables Sf ,n trade 
for all p and H. A reduction in the number of necessary variables could be achieved by introducing a default free 
referential counterparty rcf that would be the opposite side in all transactions; in this setting, the variables of interest 
would be ,rc/ , with rcf fixed, or simply S? , with an indexation of /5 such that 0 € {0, • • • ,m}. 



-^jr^ follows a beta distribution with density 

r**>M**> (*> = ^SjffgS" 1 , 0 < * < i./"^'^ («> = o */ n > i 

whore fi is the classical Beta Function 
In a preferred embodiment of the present invention, the beta family has the additional advantage of easily computed 
moments so that the moment generating function of f a%h (x) is MGf a ' b (n, a, 6) = JJ g^fcjjg and as a result calibration 
of the credit risk model is facilitated on a variety of input data types. 

s 0 

In an embodiment of the present invention, the ratios may be taken to be independent from one another, the 

*< 

beta distribution may depend upon the realized value of other inputs, where said other inputs may include different 
undcrlyings. 

In yet another embodiment of the present invention the credit risk limit of a given counterparty can inferred by setting 
a maximum responsive to the difference between the value of the counterparty liability not inclusive of credit risk and 
the value of said liability inclusive of credit risk 
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4 General Contracts Definitions 

Wc define as (/) 2 a so-called derivative security representing a contract entered into at time *o and paying 

/ (Sit * • >Sn * units °* underlying currency S c at time tn in exchange for the payment of a pr emium 

at time This more general definition includes but is not limited to the one in FAS 133 pp3-7, paragraph 6- 
11 and its 10(b) amendment in FAS 138. A derivatives contract defined with its payout payment function in the 
/ (Si , • * • , ST\ S?) format is called a Derivatives Contract Without Optional Features or DCWOF. 

Wc note: njj (/) (S\\ • • • ; S n ) = f (SiJ ■ * • \S n ) the price at time tn of the security contract entered into at time to 
and guaranteeing receipt of / (S x \ • ■ • ; S n ) 5= / (S\ 9 • - • , • ■ • , S* • • • , 5£*) units of underlying currency Sc at time tn. 
II£ ^ (/) is the trivially associated contract. 

For any i > 0 , wo define as ng >i|n (/) (S\ ; • • • ; Si) the price at time U of the security contract entered into at time to 
and guaranteeing receipt of / {S\ , • • • » S7\ • • - , fij - • , S^ 1 ) units of underlying currency Sc at time tn; ng.i.n (/) » tnc 
associated contract. 

For any i > 0 , wc define as n§ >4 (/) (Si ; * • • ; &) the price at time U of the security contract entered into at time to and 
guaranteeing receipt of ng ji+1 (/) (Si; • • • ; Si+i) units of underlying currency S c at time ng j4 (/) is the associated 
contract. 

4.1 Common Examples 

A few examples show how some well-known derivative contracts arc translated into this formal functional process: 
Vanilla Option, f (S 0 ; • - ; S n ) = (5(S n - if n )) + 

DouMe Barrier Option, / (So; • • • ; S„) = (<5(S„ - K))+ l {L <s»<H} x - - - x 1^£<si<h> 
ibian Option, / (So; • • • ; S„) = S fc ) - -FT)) 

VUatfttp St^ap, / (S 0 ; • - ; S„) = fc £ (^sfe] - I {^^^^ " 

5 Theorem of Decomposition in Basis Securities 
5.1 Derivation 

Wc have the following result further detailed in appendix and valid under the notational assumptions earlier introduced 
for any function g of m 3 : 

Theorem: For any function g defined on a discrete space of real numbers /«, with J x = I x i x • • • x J^m an m-dimcnsional 
discrete space so that I x i — {-co < x \ < • • • < x% < • • • < x l hi < +00} . With the discrete definitions of derivatives and 

2 This more general definition includes but is not limited to the one in FAS 133 pp3-7, paragraph 6-11 and its 10(b) amendment in FAS 
138. 

3 See demonstration in appendix 
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integrals earlier introduced, the following formula holds: 

5 (*V.,*") = E E E E(-i) a+p E x^x..xx*'x 

fc=0 0<< i <..<« fc <m 0<ii<-.<i p <m *=0 0<fci<..<fe,<m 

i, €{ V-.m }\{«! > fc, 6<il } 

J ( SB n{t<> - = xfiV " «))* • (*»(* - S) V -*))* 
and with Sgri{x) = 1 if x > 0 otherwise Sgn(x) = — 1 

Thus, replacing <7 by Ilg^ (/) (5b; • • • *, S») above and using our definition of basis instruments, with /? being the buyer 
and Q the seller, we deduce that, for i > 0, 

n8^i{/)(ft;.;ft_ 1 )=E E E* E * xP&J/ E E E E 

fc =00<, 1 <.<, fc <n,p=0 OSix< ; <«m A ^ O O fe= -i 



n8,*-i(/)(«i;-;*-i)-»£ E E E (-i)' +p E *? l x.x*f«x 

Jt=0 0<< 1 <.<t fc <m 0<ij<.<i p <m ff=0 0<*i<.<fcq£m 

il€{l.-»m>\{<i l ... l <j b > *|€{il,.Jp} 

i=i,.,p i*=i..,v 



(18) 



&BS ti -l ((ti,.,u),(5pn(t^ ,.,Sgn (f* - «£«■)), (F<\ ., *T), (\ju ..-Jp}\{*i. -,*•») 



Proceeding recursively backwards from n n (/) (Si; • • • ;S«) = f (S 2 ; • • • ; S n ) , we obtain IIg >0 (/) and can deduce 

o 

that II£ 0 n (/) (ft; • • • ; S n ) — E ng (i (/) (Si; • • • ; Si) The theorem of static replication of derivative securities from 

i=n— 1 

primitive securities thus follows as: 

5.2 Theorem 1 of decomposition (Static Replication) of derivatives 

Assuming no arbitrage opportunities 4 exist, any derivative security n§ >i|Tl (/) representing a contract entered into at 
time to and paying / [S\ , • • • , Sf \ - • • , S* • • • , S£*) at time *n for a premium paid at time U can be decomposed as a sum 
of basis securities as: 

j-Tl-1 

step — 1 

And wc also have: H§ A „ (/) (5 oi • ■ • ; 5 t ) = ng iS (/) (5 0 ; • • • ; 
With: 



n8j(/)(s 0 ;- -;Sj)= E . E , E E (-i) ,+p . E ^.x-xfj 



*=s0 0<»i<...« fc <m 0<ii<...<i p <m ff=0 0<fci<-..<Jb Q <m 



IbSj ((^-...^.(S^n^ 1 -^."^ il\ft.-»W)) (Soj — jS,-) 

At the end of each trading period, the sum of all debit positions is exactly equal to the sum of all credit positions. With 
positions automatically netted, this means any derivative security would be statically replicated with the selected basis 

4 The equivalence between absence of arbitrage and equivalent martingale measure is not assumed in this presentation because the possible 
existence of slippage effects and its implications for the non linearity of trading strategies break the proof of the result as detailed in [34] ,[35]. 
As a result, the whole martingale approach is not used in this presentation. 
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securities. 



The first components in formula (18) arc what we call zcroth order basis instruments "on the basis currency and arc or 
extend what is known as zero coupon bonds in a single future period case. 

Then we have zcroth order basis instruments on underlying currencies that arc or extend what is known as forwards or 
futures in a single future period case. 

First order basis instruments arc extensions on what is known as calls or puts on a single source of uncertainty in a 
single future period case. 

Second order basis instruments arc extensions on what is known as correlation options on two sources of uncertainty in 
a single future period case. 

More generally, n-th order basis instruments involve contingencies on n sources of uncertainty. 

The formula introduces those higher order terms incrementally and their premium value ends up being lower than those 
of lower order with similar contingencies. In some approximate approaches, terms of a higher order than a given value 
may be proven to be negligible. 

When the method of the present invention is merely used for pricing purpose however, it may be advantageous to use 
a basis set that makes the best use of the more directly available information. In this case the most directly available 
information is the conditional density and the associated basis set is the Arrow-Dcbrcu state-prices. 



5.3 Examples: 

1. Our World scries example is already a simple application of the definitions and formula introduced and their 
practical use outside the realm of financial assets and shows the broad scope of applicability of our method. 

2. Standard Example 

Wo will further detail how the formula works on a vanilla option example. Wo will consider the case of a single underlying 
S} = 3 t and the reference currency is iS? = 1. The formula above reduces to: 



0 

+ / 9 2 n8 j+i i Sji Fx)Bl tij , tj „ (i, -l, F j+1 x+)dx 

0 



(20) 



The derivatives of interest here is given as a vanilla call option contracted at time *o and maturing at time t n . We 
assume there is no interest rate and that the rcpo or dividend rate is zero. There is no slippage and there is no credit 
risk, but there arc n+ltrading periods *bi*i» >*n- We also assume that first order basis instruments prices arc given by 
the Black Scholcs formula 
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do = d x - crV(*i-*i-i) 

/ +»/(*>) 



With: 



<r = 0.1 
*e{l,-«,n} 
n= 1,2,7 
ti = t/365 

&toM-iM = 1 

5 0 = 100,Jr = 100 
Our goal will be to calculate ng )0 (/) as a real number. 

This example helps illustrate the impact of discrete trading against the closed form continuous trading formula. For 
the integral, we will take J x = {0, i/p, ...,p/p = 1} withp = 35, 70, 100. The results arc provided in the following 
table : 



n 


Algorithm result 


Black Scholcs 
Closed form 
continuous-time results 


1 


0.208816 


0.208816 


2 


0.29688 


0.29531 


7 


0.528928 


0.552471 


Table 2: p=35 


n 


Algorithm result 


Black Scholcs 
Closed form 
continuous-time results 


1 


0.208816 


0.208816 


2 


0.295824 


0.29531 


7 


0.533565 


0.552471 



Table 3: p=70 

These results indicate that the formula is computable and yields reasonable numbers when compared to what is known. 



5.4 Theorem 2 of compression to the formal form 

This part addresses the issue of static replication and pricing of derivatives contracts in which multiple choices arc given 
to the buyer and/or the seller of the contract throughout its life. This is a notoriously difficult problem as pointed out 
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Black Scholcs 


n 


Algorithm result 


Closed form 






continuous- time results 


1 


0.208816 


0.208816 


2 


0.295517 


0.29531 


7 


0.0.535611 


0.552471 



Table 4: p=100 



for example in US 6,321,212B1, Column 4, (7). 

The inventor of US 6,321,212B1 however docs not provide a conclusive solution, nor do other authors aware of the state 
of the art. 

Our invention, in one of its simplest embodiments teaches a solution to this acknowledged problem. Other embodiments 
solve the problem under even more complex constraints. Such constraints include accounting for microstructural issues, 
hedging strategics, etc. 

This method of our invention further enables an easily reproducible, broadly applicable and computer implcmcntablc 
sequence to solve the problem. Further still, our invention enables and teaches even more advanced hedging strategics 
and other more sophisticated systems such as exchange systems for any types of derivatives. 

Theorem 2 Let's suppose we arc in a market exchange system en tn where all possible basis instruments trade. Let's also 
suppose there arc no arbitrage opportunities. 

Lot's consider ng < n (/) a contract entered into at time to and stipulating that at time U , the seller H shall pay the 
buyer 0 f (S 0 ; A,; flo; • • • ; S n ; p n ; O,) (resp. f (S 0 \ Qq] At, ■ • - ; S n -i\ 0»; AO) 6 , where & (rcsp. H;) is an m+1- dimen- 
sional vector of parameters choices given to the buyer p at time tj and then Qj (rcsp. 0j) is the m+1- dimensional 
vector of parameters choices given to the seller ft (rcsp. /?) at time tj. Let's suppose the conditional expectations of f 
for any time in the future as a function of (Clj) 0 < j < n can be reduced to be functions with a maximum (rcsp. 
minimum) in the subset 3fcJ£ +1 (rcsp. SEn* 1 ) of possible values of ft (rcsp. Q*) 0 < j < n. 

If we have the following equality, 
/or 0 <p< iand for any function N. 

which means that if the price market makers quote for basis instruments is a detenninistic function of their parameters, 
that is the realized values of undcrlyings under consideration, 

Then ng >i>n (/) can be reduced to a contract that pays a function / (Sn; • • • ; S n ) at time £„. 
Lemma 

5 More entangled combinations of the order of choices given to the buyer and the seller can be inferred without any methodological change. 
The description outlined is merely for presentational readability. 
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If wc have the foUowing equality, 
N (Stoi • • • ftp! ftw-ii " • • » «Sti) & ^S»,ti,tj • * ' » *fc)> * * * > Wi C^i> ' ' * >-^fc)> Oit * • * i Jp)) (fto5 -ft,,; ••,5t < ) = 

/or 0 < p < i and for any function N, 
then, for any function / (S 0 \ • • • ; S n ), 

nS (i| „ (/) (ft; - • - ft; Sp+i - . . ; ft) = n£ i|tl (/ p ) (5p+i; • - • ; ft) 
forO <p< i and where f p (S p+ r, • • • ; ft) = /(fti • • - ft; ft+i; • • • ; Si) 
Proof This is a direct consequence of the decomposition formula. 

Theorem 2 proof : wc take: Ti^ = ArgMax 11% nn (— /) (ft; A; Qq\ ■ • •; ft;Ai; ^n) so that is afunctionof ft;A>; n 0 ; • ■ • ;ft»;Ai 
and /5 n = ArgMax ng fW|W (/) (ft; A)J *>o; • ■ • ; ft; A»i On) so that A* is a function of 3d; /folio;- •-;&». 
Continuing the process itcrativcly backwards, wc have for any i = n—l to 0: 

^ = ArgMax ng An (-/) (ft; A; Ho; ■ ■ ■ ; ft; A; ft+i; A+i; • • ; ft*; A>; 00 

so that J2» is a function of ft; A)> # " > fti /^»* 



A = ^irpMax ng^ n (/) (ft,; A; fi 0 ; • • • ; ft; A; ft+i; A+i ; Ei+i; ■ • * ; ft*; ftn) 

so that A is a function of ft) ; A)> * * • ; ft- 

In the end wc have the sought after function / (So; • • • ; 5 n ) as: 

/ (ft,; A,; a 0 ; • • • ; ft ; A ; ^; • • s ft; A.; ^n) 

Note: If there arc several choices for ArgMax, in any of the formulas above, any one of them can be taken randomly. 

While the ArgMin and ArgMax in the determination of Qiandfii above may not seem at first readily tractable, there 
arc numerous mathematical results that help reduce this function into boolean operators on the value of only a few 
expressions. As a result, in practice, this reverse chronological order iterative sequence may be readily implemented in 
a computerized processing system. 

Section 5.5 below further illustrates how this practical simplification is done on the most common examples. 

5.5 Examples with optimization feature 
American Options 

An American (or Bcrmudan) option here is defined here as the contract guaranteeing to the buyer the right to purchase 
(for a call) or sell (for a put) the underlying S at the strike price K at any time between t\ and t„. 

The payoff at maturity t n can therefore be written in a functional form as : 

f N x ft (tf(ft - K)) + B?, 1|n + N x (1 - A)A (*(ft - *0) + "8* Bl 2 , n + . . . ^ 
/ (ft; A; fto; ■ • ; ft; Ai On) = +ivx (l-A)-.(l- A-a)A-i (*(ft-i - *Q) + 

; +iVx(l-A)...(l-A l -0Wft-K-)) + 
with 5 = 1 for a call and 6 = —1 for a put; 0<fii<l;N being thenotionalof the contract. 
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Applying theorem 2 and given the fact that the derivative of f as a function of the is different from zero when varies 
between zero and one, we deduce that: 

A = ArgMax II§ ii>yi (/) (S Q \ A>; • ■ • ; S*; A; A+i:7*fii • • • ; S n ;P n ) 
O<0i<l _ _ 

= ArgMax IIg (i|Tl (/) (5 0 ; Ar, ■ * • I Stf A; ft+i ; A i+ i ; - • ♦ ; S«; A J 

Which can be more precisely written as: 

8=1 v ' 

^x(l~^ 1 )...(l-^_ 1 )(l~A + i)...(l-^„- 2 )/^-i(*(5„-i-K))+®BP f ,_ ltn 
+Arx(l-^ 1 )...(l-A-i)(l-A+i)-*(l^«-i)(J(5 rt -iO) + ))(5o, -A)} 

~ ^tfx(i-0i)...(i-A_i)(*<s<-JO) + ^ 

Wx(l^ 1 )...(l-^_i)(l-^ + x)...(l-/9r.-a)^i(«(5„_i-A:))+®BS^ 1 ^ 

Passport Options (One Underlying Case- Discrete translation of Hycr-Lipton-Pugachcvsky definition 6 [41]) 

A passport option is a contract giving the right to the buyer to buy or sell at anytime between to and £n~ i a maximum 
of N units of underlying S\ if at time tn the sum of those trades generate a profit the buyer keeps it, if a loss is incurred, 
the seller pays for it. 

The payoff at maturity tn can therefore be written in a functional form as: 

/ (S Q ; Aj; *V, • • • ; S n \ An; a*) = (N(fo (Si - So) + Ai (& - Si) + . . . + An-i (S n - with for 0 < i < n-1, |A»1 < 1 

As can bo seen, f as a function of any of the At 'a is twice continuously diffcrcntiablc except at most at one point. We 
thus deduce A € {—1, 0, 1} and more precisely: 

A = ^sn[n? ( <.n((A> (A - So) + . . . + A-i (A - A-i) + (Si+i - A) + A+i ($+2 - s i+1 ) + . . . + A-i (s« - - 
(A) (A - A) + . • . +A-i (A - A-t) - (s i+l - Si) + A+i {Si+2 - 5 m ) + . . . + A-i (5« - s n -i)) + ) 

(Si, 6<+i , A+i , . . . , 5„_ i , A»- 1 , S n )] 

A = 5tsn[n§,i,„((A) (Si - So) + . . . + A-i (Si - Si-i) + (S i+ i - 50 4- A+i (S <+a - Sh-i) + - . . + Ai-i (Sn - S n -i)) + - 

(A (& - So) + . . • + A-l ~ Sl-l) - (Si+l - S<) + A+l (Si+2 - S<+l) 4- . . . + A-i (Sn ~ S«-l)) + )] (&,-••. &) 

with sign(x) == 1 i/a; > 0, 5ipn(a:) = — 1 ifx < 0 . 
Passport Options (m undcrlyings case - discrete translation of Wilmott's definition 7 ) 

The passport with m undcrlyings is a contract giving the right to the buyer to buy or sell at anytime between to and 
i Tl _ 1 a maximum of N units of underlying of each of the undcrlyings S J ; if at time t n the sum of those trades generate 
a profit the buyer keeps it, if a loss is incurred, the seller pays for it. The payoff at maturity tn is: 

/(5o;A);«o;---;^n;A«;a») =JV (< Ao,(Si -So) > + <Ai,(S 2 - Si) >+...+ <An-i,(S„- Sn-i) >) + 
with 0 < i < n - 1, for Sup \p{ I < 1. 

Extending the analytical approach of the one underlying case to the multi-dimensional case, we deduce: 



A = ^rpMaxnf |4|n (/(A;A;--;A;Ai;--;^;AJ) 



6 Hyer, T. Liptom-Lifechitz,T and Pugachevsky, D. (1997) "Passport to success", Risk Magazine 10., No 9., ppl27-131. 
7 Hyungsok, A. Penaud,A. and Paul Wilxnott,P. "Various passport options and their valuation" (1999) OCIAM Oxford University Working 
paper. (42] [56] 
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" 3K£ +1 = (ZaX-.-xZar 

in which the * means at the exclusion of {0, • • • , 0} € 3fc m+1 and 2* = {-1, 0, 1} C Z 

showing that to compute , at most 3™ — 1 evaluations of expectations and comparisons will have to be made. 
Moving Window Asian option with early exercise (Hawaiian option 8 ) 

A Moving Window Asian option with early exercise is a contract in which the buyer may choose to receive at any time 
after the contract date to and the maturity tn , the difference between the average of the underlying S over the last p 
trading periods and the strike K if this difference is positive. 



The payoff at maturity tn can therefore be written in a functional form as: 

f ft (8(Si - ® Bl n + (1 - A)A - K 2 )) + <8> 5§ >n 

+ . . . + (1 - A) ■ . • U - A*-a)A*-i - €> B°-i,„ 

\ +(1 " A) • • • (1 - At-l) (*(«. - *»)) + 



/(So; A>; Ho;- •-;£«; A^OO = y x 



w£/i 5 = 1 /or a ca/Z and 5 = — 1 for a put; 0 < ft < l\p being the number of time — series 

t-1 

over which the average is taken. Ki = p x K — 
Applying the same arguments as with the American option, we deduce, 

^ " \*X<l-0l)...(l-&-l>(*(S\-*j)^^ 

+...+£ x<i-0i)...<i-&-x)<i-iM-i^^ 

1 ^x(i-^ i )...(i-^_ 1 )(«y(5 < -Jc < )) +v '& l, B^. i „>ng fifrt (f x(i-^i)...(x-^-i)A + i(tf(s <4 -i-^+i)) + ®BS,^i. rt +- 

+^)<(l-^ 1 )...(l-^„ 1 )(l-^, +1 )...(l-^> 2 )^„_a(tf(&' rt ^ 1 -K„„ 1 )) + ®Bg t ^ l(n 
+&x(i-A)...(i-A.i)(i-A + i)»(i-^„_i)(«(s„-ic„))+))(So t ».,S 4 )> 

6 Applications 

The methods of the present invention will be further illustrated in preferred embodiments in conjunction with three sys- 
tems of practical use for trading purposes: pricing systems, derivatives exchange systems and risk management systems. 



6.1 Derivatives Pricing 

In pricing or decision systems, the methods of the present invention will be used as alternatives to known methods with 
four main practical advantages over closed form, trees, PDB, Monte Carlo: 



• Flexibility in the design of any conceivable derivatives contract, even those where subsequent choices may be made 
by parties in the contract, that will be seamlessly priced for any specification of the dynamics of the undcrlyings. 

• Easy incorporation of liquidity effects in the pricing of derivatives through the introduction of scaling density 
functions. 

• Easy incorporation of credit risk through tho seamless integration of undcrlyings representative of the counterparty 
credit risk. 

8 wihnott.com, Technical forum (2001) 



WO 03/107137 



37 



PCT/US03/19179 



• Deriving hedging strategics independent of model risk through the use of indicated basis instruments contracts 
instead of Greeks as is the current practice derived from Black-Scholcs- Mcrton. 

As far as speed is concerned, except for closed form solutions, the methods of the present invention compare to those 
of PDEs, Trees or MontcCarlo and many of the speed improvement methods of regular integral may be advantageously 
integrated in this framework; such methods would include Simpson types of methods, Gaussian quadratures and related 
methods, importance sampling methods, low discrepancy methods on Sobol, Halton or other family of points. 

Sec further details in the drawings. 

6.2 Exchange Design 

Many types of derivatives exchanges may be derived from the methods of the present invention. We present here an 
example of design that seeks to extend the current practices and derivatives contracts to achieve the uses outlined for 
the present invention. However the example can be readily applied in any other OTC or regulated exchange framework 

Sec further details in the drawings. 

6.3 Risk Management Systems 

In derivatives portfolios Risk Management, traders arc given Value At Risk (VAR) and/or Greek limits to reduce port- 
folios exposure, to adverse market conditions. Distributional assumptions used to infer VAR values as well polynomial 
approximations implied by Greeks use create model risks that arc poorly understood in this framework. 

The existence of a market for basis instruments implies assumptions-free market state price densities for which instru- 
ments exist to manage any risk ono would want to hedge. 

Going beyond Value At Risk, in an exchange market of basis instruments, more specifically tailored derivatives may 
help better manage portfolios of derivatives. 

The methods of the present invention can be used to design and build more effective risk management systems. In a 
prcforrcd embodiment of the present invention, such risk management systems include: 

• Databases of derivatives contracts positions stored in basis instruments units 

• Inputs feeds of market or simulated basis instruments prices 

Such a system may then be connected to an exchange system as specified above and enable effective trading strategics 
to match the desired outcome of the risk manager. 

For instance a portfolio manager with a portfolio of assets including assets derivatives with pay-offs at fk(S 0 ; S\\ • • • ; S nk ) 
for fc = 0, ,p and wishing be hedged in a value at risk framework by purchasing the contract ng >i|<+J .( £ f k ,VAR(i,i+j)) 
entered into at time *o hi which the assets portfolio manager would bo paying at time U the premium 

nS A4+J ((nS,i(£ /fc(5 0 ;S i; • - , ; S nh )) - IIS, i+i (£ f k (S 0 ; ft; • • • ; S nh ) + VAR(i,i + j))+) 

fc=0 k=0 
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and receiving at time the quantity 

. So if wc have 1 — cvAnfaj) the confidence level associated with the amount VARfai + j) , then cvAR&j) is 
characterized by: 

yAR(i,i + j)=Min{5>0.sticAt/wii(£7(l r • ^ o o xxm ,J> e e 0 ^ v )) < cvAii(*,i)} 
v " 1 WjIE MSoiSn ••;5„ fc ))-n« > ^(E /*(So;Sii -;S» fc )>-*} 

A.— O fc=0 

The main drawback of the value at risk method that is resolved in another preferred embodiment of the present invention 
is the possibility to concentrate or hide excessively negative payoffs in states where the potential loss exceeds the com- 
puted value at risk amount. Fig 7 shows a better illustration of this problem. In a risk management approach consistent 
with methods of the present invention, a custom portfolio payoff may be designed to provide a maximum possible loss 
for all possible cases. The method may them be implemented in risk management systems or risk management computer 
program products. 

Sec further details in the drawings. 

The method may them be implemented in risk management systems or risk management computer program products. 

In yet another embodiment of the present invention derivatives accounting responsive to FAS 133 or IAS 39 is facilitated 
by removing uncertainty that comes with assuming a particular model to allocate profit or losses on derivatives as 
hedging instrument and hence their allocation an Other Comprehensive Income category. As such it is more effective 
than the methods disclosed in WO 02/44847 A2: Dynamic reallocation hedge accounting, 
US20020111891 : Accounting system for dynamic state of the portfolio reporting, 

US20020 107774 : Compensatory ratio hedging. All those methods arc hereby incorporated by reference. 
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7 Dimension & Computations Reduction through Changes of Variables 

In a practical implementation, the computational time is a major concern. It may appear at first computationally in- 
tractable to implement a full-dcpcndancc for the price of the option of all the So, . . . , Sn because otherwise the algorithm 
computation time grows exponentially with the number of time increments. 

In practice, however, for virtualy all known payofls of interest or basis instrument contracts used, relatively obvious 
changes of variable may be made that result in the computational time ending up in the order of a power two or three 
of the number of time steps. We review some of them below. 

The purpose of the algorithm is so to decrease the complexity of the option price dependence. It can be easily ac- 
complished for non-path-dependant options and undcrlyings following a martingale, hence markovian process ( process 
driven by brownian motions or more generally Levy processes). For example, the price of an European Vanilla (under 
geometric brownian motion assumptions) at time t only depends of the curcnt value of the underlying St and not on its 
past values. In these cases, we consider the option price as a fonction of the only current underlying value. For example, 
we have : 

n§, i( „(/)(5o, • • • , = ns ii>n (/)(£) 

With such an implementation, the computation time is almost linear in maturity 9 . But we can't apply this method to 
more general payoff that arc path-dependant such as asian options, loop-back options, amcrican options, etc... 

To implement this type of path-dependant option, we consider the price of the option as a function of several variables 
V?'. one will be the current value of the underlying and the other variables will be used to compress the path-information 
needed to determine the price of the option at any time. For example, we choose for: 

• the Arithmetic asian option: we use a two-dimension variable with = Si and = Z)]Uo ^fc- We define 
the payoff function as a function of the vector: f n (V) = (V£ — iif) + . 

• the loopback option: we use a two-dimension variable with V? — Si and V^ 1 = max^ =0 S*. We define the payoff 
function as a function of the vector: f n (V) = {V£ — K) + . 

• the Amcrican option: we use a two-dimension variable with V? = 5» and V* is used to store the information that 
the option has been or not exercised in the past and in the first case the profit that have been realised. That is: 

yi _ \ 0 if the option is yet to be exercised 

IP > 0 if the option as been exercised with a profit of $ p 

We define the payoff function as a function of the vector: f n (V) = (V£ — Jf) + . 

• the Passport Option: we use a two-dimension variable with V? = Si and is used to store the current profit of 
the process. That is: 

We define the payoff function as a function of the vector: f n {V) = V^J. 

With this method, we need to implement the link between these variable, it means that, given the vector of variables at 
time U and the value of the underlying at time U+i* how to determine the value of the vector of variables at time 
In the past examples, this relation is obtained as follow: 

9 The implementation is actually notlinear in time because the range of points to consider increases with time. 
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the Arithmetic asian option: 



• the loopback option: 



1 i+1 i + i 1 



the American option: at time t £+ i, we have two possibihtics, cither we have already exercised the option so nothing 
to be done, or the option is yet to exercised and we need to make the choice of exercising it. We decide that the 
option will be exercised if the resulting profit is higher than the future expectation of profit at maturity of this 
option that is the price of the European option with spot value the current value of the underlying with strike and 
time of maturity similar to the American option. That is, for the call particularly: 



V> if V; 1 > 0 



if V? = 0 then 



Calln(Si+u K) - (S i+ i - K)+ if 
0 if 



CalUSi+u K) - -K)+ > 0 
Call n (S i+u K) - (S i+ i - K)+ < 0 



For the practical implementation, wo first compute with our algorithm the price of the corresponding European 
option that gives the price of the European option at any time for any spot equal to the past values of the 
underlying which is exactly what wc need to compute the algorithm for the American Option. 



r 


b 


a 


n 


bp 


prcc 


Result 


Bjerksund and Stcnsland 10 
continuous- time 
Approximation 


European 
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-0.1 


0.1 
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0.05 
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0.4757 


0.4680 
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0.05 


10 


0.05 
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0.3981 


0.3974 


0.3981 


0.1 
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0.1 


10 


0.1 
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0.7953 


0.7918 


0.7950 


0.025 


-0.05 


0.05 
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0.1 


5 


0.3578 


0.3159 


0.3084 
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-0.1 


0.1 


10 


0.05 


4 


0.7573 


0.7396 


0.7025 


0.05 


0.05 


0.1 


10 


0.05 


5 


0.8933 


0.8969 


0.8983 


0.05 


-0.05 


0.1 


10 


0.05 


5 


0.7357 


0.7046 


0.6989 



Table 5: American Option of spot So = 100 and strike K = 100 Maturity n in days, Interest rate (r), Cost of carry (b), 
Basis Point & Precision. 



• the Passport option: at time U+i , wc choose ft+i such as the profit A+i - Si) is maximum. The link between 
the vector at time U+i and time t» is exactly the recurrent definition of V^. v 

General Computational gains approach: 

In addition to the methods described above, in general, computational gains may be achieved by the following process: 



1. define / in discrete space 
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2. extend / in the continuous space with an infinitely smooth or analytic function fc that interpolate the values of / 
in the discrete space 

3. use well known theorems to obtain the value of the integral by using an approximating function of / and use a 
smaller subset of the original discrete space to obtain the approximate value of the integral 

This approach eliminates the need to verify the functional regularity conditions that one considers when taking from 
the outset a continuous space. It results in methods, systems and computer program products that run faster because 
larger numbers of states can be spanned at once. 
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APPENDIX 

Theorem 

For any function g defined, on a discrete space of real numbers J s , with I x = J»i x • • - x I x ™ an m-dimcnsional discrete 
space so that J x , = {-co < x{. < • - - < x 0 < - • < xj^ < +00} . With the discrete definitions of derivatives and integrals 
earlier introduced, the fbUowing formula holds: 

g(x\.. t aT) = £ E E E(-l)' 4 * E a^x..x^x 

/ (*» (<•» -*SV - •• - *)V - *)) + 

and with Sgn(x) = 1 if x > 0 otherwise Sgn(x) = —1. 

We establish the proof in a discrete setting in order to avoid making regularity assumptions that may end up not 
being true for cases of interest and to match real world discreteness. These cases of interest include the most common 
derivatives pay-off structures (S T - i0 + that is a function of St that is not twice differentiate in that variable. 
Proof: 

Case 1: One Variable 

Let's consider l^h integers such that I < 0 < h and xi,,x 0 ,,xh is an increasing sequence of /i— I + 1 real numbers. 
We note J x = {-co < x u , x 0 , ,x/ t < +00} and let's take x n e I* . For any function g defined on J x> 
If x n > xo, We have: j 

g(x n ) = gfro) + te(x n ) - g(x 0 )) 

n-1 

g(*n) = gfro) + E (5(^+1) - 9(*j)) 

i=o 

n-1 

$(*n) = 9(z 0 ) + E (*i+i " 



Likewise, we have: 



i-i 

^'(Xj) = p'(xo) + X) - *kW(**) 



So replacing in $(x n ) gives: 

y(*n) = ff(xo) + if (x i+ i - xj) (g f {x Q ) + E (**+i - afcV'O&fc) J 
i=o \ fc=o ^ / 

= g (* 0 ) + E te+i - *i)^(«o) + E - E (^fc+i - *fc)*/"(*fc) 

j=n jf=0 fc=0 

Xhcn, inverting the order of the summation on the indices k and j gives: 

g(* n ) = g(x Q ) + (x n - x 0 )g , (x Q ) + n ff "ff (x j+1 - x^(x M - x k )^{x k ) (*) 

fc=0j=fc+l 
n-2 

$(*n) = S(*q) 4" (X„ - XoV(x 0 ) 4" E (*n - *Jb+l)^"(*Jb)(*fc+l " **) 

fc=0 
n-1 

$(*n) = g(x 0 ) + (Xr, - X 0 V(X 0 ) + E (*n - Xfc+iV (xjb)(x fc+ l - Xfc) 

fc=0 

thus we can rewrite: 
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+00 



If x^ < xq, we have: 



Likewise, we have: 



So replacing in g(x n ) gives: 



9(* n ) = 0(*o) + E - sfri+i)) 

= <7(*o) + E - 

-1 



g(*n) = *(*o) + E (*i - (^(aso) + E («* -*fc+iV(*0 ] 

<K*~) = *(*o) + E (*j - sh-i)^ M + E (*i - E (*fc - 

j=t» i=n fc=i 

-1 -1 
= 9(*o) + (a^ - ar 0 )s'(*o) 4- E (*i - E (** - *fc+i)y"(a:fc) 

Then, inverting the order of the summation on the indices k and j gives: 

g(x n ) = g(x 0 ) + (x n - x 0 )g , (x 0 ) + E (*i - E (*fc - *k+i)g"(*k) (*) 

i=n fc=j 
-1 fc 
y(x n ) = p(x 0 ) + («» - SoM^o) + E (** ~ Sk+i)s"(*Jb) E (*j " 

-1 

y (x n ) = 0(x o ) + (x n - x 0 )^(a;o) + E (** - ifc+iV'Wfe - 

A: =71 

So that we can write: 

ff(*n) = tf(*o) + (a?« - aroVfco) + E g"{*k)(xk+i - x n )(x k +i - or : 

fc=n 

= 9&o) + (*n ~ 2:0)^(^0) + E - * n ) + (*fc+i - x fc ), which is also : 

= *(*o) + (*« - *oV (*o) + / 0" (*)(*+ - 

—00 

Adding the cases when x n > x 0 and x n < xo, gives: 
Lemma 

For any x e I* = {—00 < xi, ,x 0 , ,x/ t < -f 00}, with our discrete definition of derivatives and integrals, we have: 

h-2 

g(x) = ^(x 0 ) + (x - xoV(xo) + E g"(*k)(Sgn(x k - x 0 )(x - x fc+ i)) + (xjfe + i - x*), which is also : 

fc==i 

= p(x 0 ) + (x - xoV(x 0 ) + rg"(t)(Sgn(t - x 0 )(x - 

—00 

Case 2: m variables m > 1 
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a) Preliminaries 

First wo introduce objects and definitions customized from the theory of distributions. Wc consider the m-dimcnsional 
discrete space I w = I x i x • • • x l x m , with I x i = {-co < a;} < • • • < xj < • • • < x^ < +00} 

S x t is an application such that for any application g defined on a subspacc of I x , 

S xh : g {x\-.-,x i -\xW+\---,x m )~*g(x\---,x i -\xi ) ,x i +\---,x m ) = 
< S a:i) ,g(x 1 ,- -,x i -\x i ,x i+1 , --,x m ) > 

S k -9( xl > ~> **>* <+1 . — ,*")-* -dig (x 1 , ■ ■ ■ , a*~ l , xj, x* 1 , • • - ', x m ) 
=< , g(x\..., x'-'.x', x i+ \ • • • , x m ) > 

C a ■■9(x\---,x i - 1 ,x i ,x i +\..-,x m )^dfg(x l ,- - t x i -\x} ) , x™ , • • • , x"») 

Wc also introduce the operator of composition o of the£ s <, , (5^/s as: 

S <4 oS *o ■■9(*\---,* i -\xW + \- -,J-\*,^ + \---,x m )~g(*\---, x 0 , x** 1 , • • • , xS-\4, x>+\ • • • , x™) 
=< S^oS^g (*\ • • • , x*- 1 , x*. x* +1 , • • • , a*" 1 , x>\ x^ 1 , ■•■,x m ) 

5 k° S 4 

= S^oS'^ : g (x l , • • • , x*" 1 , x\ x i+1 , • - • , »*-», x>, • • • , x*») - (x 1 , • • - , x*-\ x«, x<+\ • • • , x*" 1 , xj, x*+\ • • • , x m ) 
=<*-i otf ii.»( !cl . ••,x < - 1 1 a i ,x < + 1 ,-..,x»- 1 ,x»',xJ +1 ,.-. 1 x TO ) > 

= ^< 4 :ff(x 1 ,-", a : < - 1 ,x*,x*+J, ...,xJ-^x^,xi+^...,x™)^a? ^ ,(x^.••,x*-^ a 4,x^^...,xJ-^xS,x^+^...,x™) 
=<^oo^,»(* 1 ,--,a: 1 - 1 ,x < ,x*+S...,xJ-\xi,x*+i..-. | x»*)> 

So, g (x\ • • • ,x m ) =< o^n-o • • . o5 xl ,g (t 1 , •••,*"*)>. 
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From these notations, it becomes obvious that the result obtained in our lemma can be rewritten as: 

< <5x, 9(v) >=< 4e 0 . > -(» - *o) < S '*o • > + / ( 5 *»( x - " x )> + < 5 « » ff(u) > di 

J* 

or rrwre conveniently, 6* = £ xo - (a: - ar 0 )<% 0 + / (Sgn(x - a?o)(*+ ~ 
b) Proof by recurrence 

Since it clearly appears that this operator is distributive with respect to the addition, we will assume that for m£l, 
the following relationship holds: 

5^- l0 ...o«,x= E 1 E E feW* . 

JkssO 0<< 1 <...<t ij <m-i 0<ii<-.<ip<m-l g«0 0<fci<...<fc 9 <m-l 

Ji€ii.".m-i}\{<i fc|€iil..--Jp> 

1=1, ...,p 1=1, ...s 

dt il -di ih with {*i,...,*r} = {1 m-l}\-(ti > ...,tib > ii l ...Jp} 

Now we must prove the formula for £r«»o • • - o*«i . From the one variable case we know that: 

6 m m = S x r-(x m -aK)8' xh n+ j (Sgn{t m -xZ)(x m -tZ)) + 6" m dt m 

So, 

S^oS^-io^oS^ = ^ - (x m - xZK~ + / (Sffn(i m ~aff)(a: m -^^ 

^m-io-^^E 1 E E £(-1)^ E ^x-.x^-x 

k=0 0£4 1 <...«*£m-l 0<i 1 <.-.<i p Sm-l 0=0 0<fej <...<fc,<m-l 

Jl€{l.-,m-l}U*i <fc> *i6{ii,.. ; Jp} 

i»i,...,p i=i,--»« 
o - (x m - rrS 1 )^ + / (Sgn(t m - xf)(x m - i7)) + <£»<ft m J dt'* - • • <&» with {si *•} = {!, .... m - , t* ( 
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fc=0 0<<i<...«j,<m-l 0<ii<...<i P <m-l 9=0 0<kK ; ..<k,Sm-l 

l=l»— .p 1=1.— .9 

/ (5 fl n(^-4 l )( a:<l -*i l )) + -( 5 fl"( t4fc -^) (* <k -*+)) + fe^ x^°" - -<Jp^° 

oi.g.dt* 4 •••di lfc + 

m t l E E f:(-D 9+ ^ 1 E xj'x-x^x 

Jfe=0 0<u<...<i fc S»»-l 0<ii<...<i P <m-I 9=n 0£fci<...<fcg^m-l 

/ (** ~ *m**- • • • fa - 4 fc ) («* : *J )) + ^r^:x m * 

5^ o • • • o^- fc o*^ o • • • o&j, o^p oS^ <> • • o5 x . r dt*i • - ♦ dt <fc + 

m f? £ £ £(-l) ff+p E ^x-x^x^x 

Jk=0 0<<i<...<t A .<m-l 0<iK...<ip<m-l 9=0 0<ki< ...<fc*£m-l 

i|€{l."<m-l}\{<i fcie«i,...Jp) 

f=l.- ,p I=l,—,9 ^ 

/ (Spn - 40 (* <l " • " • (t'» - 4 fc ) (** " 

m E E E E(-D« +p E -j»x...x«frx 

Jb=0 0<<i<...<t fc <m-l 0<ii<...<i p <m-l C=<> 0<fci<...<fc,^m-l 

iie{l,-.m-l>Mii.....*Jb} fc|€-lJi,...JpV 

(00» (t*» - 40 - iV')) + • • • (S^n (t<> - 4") (* <fc - #)) + (5 9 n(£ m - affX*" - t?)) + feffifr * 
5" 4l o- • • od'^aSlmoi^o • • • oS'^oS^o ■ ■ ■ oS x .-r eft* 1 • • dt<*dt m with {s L , s r } = {1, ...,ro - l}\{ii, ...,<*, Ji, .... i P } 

(-0 = "E E E E(-1)' +p E ■8»X — x«S«X 

fcz=0 0<ti<...<tfc^m-l 0<i 1 <...<ip^m-l 9=0 0<fci<...<k q £m-l 

/ (5pn (t'i - 40 («<• - 't{)) + ■■ ■ (Sgn (*> - 4") - fe*f£ x " 

di<i . . . dt 1 * tirit/i {si, ...,s,} = ({1, ...,m - -,**. ill --.ip}) U {m} 



(") = "E E E £(-D' +,> E *8»x...x«frx 

A=0 0«x<...« fc <m-l 0<ii<...<ip<m 9=0 0£* 1 <...<* y <m 

Jl€«Y".m-a>M<i.... l *Jb>)tlim> ^^^x^a 

dt fl • . with {s u «r } = ({1, -.,971- Jli "iip» 



V E E E(-D' +p E o4»x.-xa^x 

fc=0 0£<x<...«hSm-l 0<ii<-..<i p <m-i 9=0 0£fci<...<fe g ^m 

/ (t*i - 40 " - (Sgn (#■ - 4") (x<- - $)f x 

dt* 1 tutf/i {*i,...,sr} = ({l,.. M m- l}\{t*i ii» -iip}) 
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W=E l E E E(-i) e+p E ^'x.-x^'x 

/t=Q 0«i<...<t fc «=m 0<iK..-<ip<m-l ff~0 0£k 1 <...<* g <m-l 

i!6<l,-.m-l>\{<i,...»*fc> fc|€{ii,...Jpf 

Thus, summing (I), (II), (III), and (IV) , we can easily see that can be rewritten as: 



Jb=K) 0<ti<...«fc<m 0<i 1 <...<i p Sm 9=0 0<*i <...<k 9 <m 

il€{l, -.m}\iii i k } fcie{ii,-..Jp> 

• / (Spn (a * -^)) + ...(5sn(t"> -afr) (*<> -#)) + fs£i£|£ X ^o- o^o^o- • 

dt il - dt ih with {ai,..., s r } = {I, ... ,m}\{ti,. t fc| 

which is the result sought to prove our formula by recurrence. 

In a continuous space representation and if / is twice derivable, m = 1, t+ is replaced by t and we recover the well 
known formula derived using other tools for instance in [13]. 
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All patcn&citcd in this application arc hereby included by reference. 

While the present invention has been described in connection with preferred embodiments, it is not intended to limit 
the scope of the invention to the particular form set forth, but, on the contrary, it is intended to cover such alternatives, 
modifications, equivalents as may be included within the spirit and scope of the invention as denned in the appended 
claims. 



